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ABSTRACT 

Context. Chamaeleon I is the most active region in terms of star formation in the Chamaeleon molecular cloud complex. Although it 
is one of the nearest low-mass star forming regions, its population of prestellar and protostellar cores is not known and a controversy 
exists concerning its history of star formation. 

Aims. Our goal is to search for prestellar and protostellar cores and characterize the earliest stages of star formation in this cloud. 
Methods. We used the bolometer array LABOCA at the APEX telescope to map the cloud in dust continuum emission at 870 yum with 
a high sensitivity. This deep, unbiased survey was performed based on an extinction map derived from 2MASS data. The 870 yum map 
is compared with the extinction map and C'^'O observations, and decomposed with a multiresolution algorithm. The extracted sources 
are analysed by carefully taking into account the spatial filtering inherent in the data reduction process. A search for associations with 
young stellar objects is performed using Spitzer data and the SIMBAD database. 

Results. Most of the detected 870 fim emission is distributed in five filaments. We identify 59 starless cores, one candidate first 
hydrostatic core, and 21 sources associated with more evolved young stellar objects. The estimated 90% completeness limit of our 
survey is 0.22 Mq for the starless cores. The latter are only found above a visual extinction threshold of 5 mag. They are less dense 
than those detected in other nearby molecular clouds by a factor of a few on average, maybe because of the better sensitivity of 
our survey. The core mass distribution is consistent with the IMF at the high-mass end but is overpopulated at the low-mass end. In 
addition, at most 17% of the cores have a mass larger than the critical Bonnor-Ebert mass. Both results suggest that a large fraction of 
the starless cores may not be prestellar in nature. Based on the census of prestellar cores. Class protostars, and more evolved young 
stellar objects, we conclude that the star formation rate has decreased with time in this cloud. 

Conclusions. The low fraction of candidate prestellar cores among the population of starless cores, the small number of Class 
protostars, the high global star formation efficiency, the decrease of the star formation rate with time, and the low mass per unit length 
of the detected filaments all suggest that we may be witnessing the end of the star formation process in Chamaeleon I. 

Key words. Stars: formation - ISM: individual objects: Chamaeleon I - ISM: structure - ISM: evolution - ISM: dust, extinction - 
Stars: protostars 



1. Introduction 

Large-scale, unbiased surveys in dust continuum emission have 
greatly improved our knowledge of the process of star forma- 
tion in nearby, star-forming molecular clouds. During the past 
decade, (sub)mm dust continuum surveys established a close 
relationship between the prestellar core mass function (CMF) 



* Based on observations carried out with the Atacama Pathfinder 
Experiment telescope (APEX). APEX is a collaboration between 
the Max-Planck Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 
** Appendices A, B, and C are only available in electronic form via 
http://www.edpsciences.org 

*** The image shown in Fig.[2]and the content of Table[2]are available 
in electronic form (FITS format for the former) at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5). 



and the stellar initial mass function (IMF), suggesting that the 
IMF is already set at t he early stage o f cloud fragmentation 
(e.g. Motte et al. 1998; iTesti & Sargent [l998; Johnstone et al] 
200(J). In these submm surveys, prestellar cores are found only 
above a visual extinction of 5-7 mag, suggesting the existence 
of an extinction threshold for star format i on to occur (see, e.g., 
Johnstone et al. 2004; Enoch et al."2006^, 'Kirk et al.' 2006, but 
also Onishi et al. 1998 based on molecular line observations). 

In the mid-infrared range, five of these nearby molecu- 
lar clouds - Chamaeleon II, Lupus, Perseus, Serpens, and 
Ophiuchus - were mapped with good angular resolution and 
high sensitivity with the Spitzer satellite in the frame of the c2d 
legacy project (From molecular cores to planet forming disks, 
see Evans et al. 2009). This extensive project produced a deep 
census of young stellar objects. It delivered sta r form ation effi- 
ciencies ranging from 3 to 6% and lEvans et al.l (l2009t) estimate 
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that these efficiencies could reach 15-30% if star formation con- 
tinues at current rates in these clouds. Star formation is found to 
be slow compared to the free -fall ti me a nd hig hly concentrated 
to regions of high extinction dEvans et alJl2009l) . consistent with 
the threshold mentioned above. 

Recently, the far-infrared/submm Herschel Space 
Observatory revealed an impressive network of parsec- 
scale filamentary structures in near by molecular clouds in the 
frame of the Gould Belt Survey dAndre et alj |5010a). This 
survey greatly increases the statistics of prestellar cores and 
confirms the close relationship between the CMF and the IMF. 
In the Aquila molecular cloud, sever al hundreds of prestellar 
cores are identified, and lAndre et"aD (l2010a) show that most 
of these prestellar cores are located in supercritical filaments. 
They confirm the existence of an extinction threshold for star 
formation at Ay ~ 1 mag and propose that it may naturally 
come from the conditions required for a filament to become 
gravitationally unstable, collapse, and fragment into prestellar 
cores (see also lAndre et al.ll2010bh . This finding will be further 
investigated in the frame of this Herschel survey, but it can also 
be partly tested in other nearby clouds in the submm range, 
thanks to the advent of very sensitive cameras like LABOCA at 
the Atacama Pathfinder Experiment telescope (APEX). 

Chamaeleon I (Cha I) is one of the nearest, l ow-mass star 
form ing regions in the southern sky (150-160 pc. Whitt et et alJ 
[1997; Knude & H0g 1998; see Appendix |RT] for more details). 
Together with Cha II and III, it belongs to the Chamaeleon cloud 
comple x whose p opulation of T-Tauri stars has been well stud- 
ied (see iLuhmanl bOOS and references therein). Cha I contains 
nearly one order of magnitude more young stars than Cha II, 
while Cha III does not have any. Surveys performed in CO and 
its isotopologues showed that Cha I contains a larger fraction of 
dense gas tha n Cha II and III alt hough the latter are somewhat 
more massive (iBoulanger etal.ll9 98; Mizuno et al. 1999, 200l|). 
Finall y, several in dications of jets an d outflows we re found in 
Cha I j Mattila et al. 1989; Gomez et al..2004; Wang & HenningI 
l2006t lBelloche et al. 2006). Only three Herbig-Haro objects are 
know n in Cha II and none has been found in Cha III ('S chwartj 
[1973). Cha I is therefore much more actively forming stars than 
Cha II and III. 

A controversy exists, however, concerning the history of star 
formation in Cha I. Based on the H-R diagram of the young 
stellar population known at that time (47 members with known 
spectral type), Palla & Stabler (2000) conclude that star forma- 
tion in Cha I began within the last 7 Myr and that its rate 
steadily increased until recently. They find similar results in sev- 
eral other clouds and draw the general conclusion that star for- 
mation in nearby associations and clusters started slowly some 
10 Myr ago and accelerated until the present epoch. They spec- 
ulate that this acceleration arises from contraction of the parent 
molecular cloud and results from star formation being a critical 
phenomenon occuring above some threshol d density. Howe ver, 
based on a larger sample of known members, Luhman (2007*) re- 
cently came to the conclusion that star formation began 3-6 Myr 
ago in Cha I and has continued to the present time at a declining 
rate. Since little is known about the earliest stages of star for- 
mation in Cha I, investigating whether there is presently a pop- 
ulation of condensations that are prestellar, i.e. bound to form 
stars, should provide strong constraints on these two competing 
scenarii for the star formation history in Cha I. 

To unveil the present status of the earliest stages of star for- 
mation in Cha I, we carried out a deep, unbiased dust contin- 
uum survey of this cloud at 870 //m with the bolometer array 
LABOCA at APEX. The observations and data reduction are de- 



scribed in Sects. |2]and[3] respectively. The maps are presented in 
Sect.m SectionOexplains the source extraction and the sources 
are analysed in Sect. [6] The implications are discussed in Sect.|7] 
Section[8]gives a summary of our results and conclusions. 

2. Observations 

2.1. 870 fim continuum observations witii APEX 

The region of Cha I with a visual extinction higher than 3 mag 
was selected on the basis of an extinction map derived from 
2MASSE1 (see Sect. It was divided into four contiguous 
fields labeled Cha-North, Cha-Center, Cha-South, and Cha- West 
with a total angular area of 1.6 deg^ (see Fig. [TJ- The four 
fields were mapped in continuum emission with the Large APEX 
Bolometer CAmera (LABOCA, ISiringo etalJ l2009) operating 
with about 250 working pixels in the 870 fim atmosphe ric win- 
dow at the APEX 12 m submillimeter telescope (Giiste lTet al.l 
2006). The central frequency of LABOCA is 345 GHz and its 
angular resolution is 19.2" (HPBW). The observations were car- 
ried out for a total of 43 hours in August, October, November, 
and December 2007 (fields Cha-Center and Cha-South) and for 
41 hours in May 2008 (fields Cha-North and Cha- West), un- 
der excellent (t^™^" = 0.13) to good (t^™^" = ^.37) atmo- 
spheric conditions. The sky opacity was measured every 1 to 
2 hours with skydips. The focus was optimised on Mars, Saturn, 
or Venus at least once per day/night. The pointing of the tele- 
scope was checked every hour on the nearby quasar PKS 1057-79 
and was found to be accurate within 1.9" for the 2007 data and 
2.8" for the 2008 data (rms). The calibration was performed with 
the secondary calibrators IRAS 13134-6264, IRAS 16293-2422, 
V883 Ori, or NG C 2071 that were o bserved every 1 to 2 hours 
(see Table A.l of lSiringo et al.ll2009l) . Measurements on the pri- 
mary calibrator Mars were also used. 

The observations were performed on the fly, alternately with 
rectangular ("OTF") and spiral scanning patterns (see Sect. 8 of 
ISiringo et al.1l2009l) . The OTF maps were alternately scanned in 
right ascension and declination, with a random position angle 
between -12° and +12° to improve the sampling and reduce 
striping eff'ects. The "spiral" maps were obtained using rasters 
of spirals (see Sect. l3.1l for the parameters), with a random offset 
around the field center (within ±30" in each direction) to also 
improve the sampling. 



2.2. Extinction map from 2MASS 

We derived an extinction map toward Cha I from the publicly 
available 2MASS point source catalog (see Fig.[T]). We selected 
all sources detected in either both J and H, or both H and K, or in 
all three bands. Before computing the extinction map, 223 young 
stellar obj ects (YSOs) th at are known members of Cha I (see 
Table 1 of lLuhmanll2008h were removed from the sample. The 
extinction was calculated from the average reddening of stars in- 
side the elements of resolu tion of the map using the method de- 
scribed in Schneid er et alJfcO lO) that is adapted from Lada et al.l 
(119941) . ILombardi & AlvesI fcoOl.) . and ICambresv et alj (|2002|) 
The extinction toward each star is obtained from the uncertainty- 
weighted combination of the [J-H] and [H-K] colors. The as- 
sumed average intrinsic colors are [J-H]o - 0.45 + 0.15, and 



The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the 
National Aeronautics and Space Administration and the National 
Science Foundation. 
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Fig. 1. Extinction map of Cha I derived from 2MASS in radio 
projection. The projection center is at (a, 5)j2(K)0 - (ll''01™24^ 
-77°15'00")- The contours start at Ay = 3 mag and increase by 
step of 3 mag. The dotted Hnes are Hnes of constant right ascen- 
sion. The angular resolution of the map (HPBW - 3') is shown 
in the upper right corner The four fields selected for mapping 
with LABOCA are delimited with dashed lines and their center 
is marked with a black cross. The white plus symbols in fields 
Cha-Center and Cha-North mark the positions of the dense cores 
Cha-MMS 1 and Cha-MMS2, respectively, detected in dust con- 
tinuum emission at 1.3 mm with the SEST (Reipurth et al. 199^. 
The field of view of LABOCA is displayed in the lower right 
corner 



[H-K]o = 0.12 + 0.05. These are derived from the typical disper- 
sions for a population of galactic stars as measured using simu- 
lations with the Besangon stellar population model (Robin et al.' 
[20030. The infrared color excess was converted to visual extinc- 
tion assuming the standard extinction law of fRieke & Lebofskvl 
with a ratio of total to selective extinction Ry - 3.1 (see 
Appendix lB.3l l. In addition, we use the Besanfon galactic model 
to derive the predicted density of foreground stars in the 2MASS 
bands. For the distance and the direction of the Chameleon 
clouds, it is found to be less than 0.01 star per arcmin-. Therefore 
it is negligible and no filtering for foreground stars is applied. 

A Gaussian weight function for the local average of the in- 
dividual extinctions defines the resolution of the final map, with 
a FWHM of 3' and a pixel size of 1.5'. The FWHM is cho- 
sen so that most pixels (97%) with Ay < 6 mag have at least 
10 stars within a radius equal to FWHM/2. Most pixels (96%) 
between 6 and 9 mag and between 9 and 13 mag have at least 
4 and 3 stars, respectively. All pixels above 15 mag have be- 
tween 2 and 4 stars. Only five pixels have no background star 
within FWHM/2. Their extinction value is solely estimated 
from the stars located between FWHM/2 and the truncation ra- 
dius (1 .3 FWHM), which means a small local loss of resolution. 



^ See also |http://www.obs-besancon.fr/| 



Two of these pixels are located within the Ay - 12 mag contour 
of the Cha-MMS 1 region. 

The typical rms noise level in the outer parts of the map is 
0.4 mag, i.e. a 3cr detection level of 1.2 mag for an FWHM of 
3'. This rms noise level is however expected to increase towards 
the higher-extinction regions because of the decreasing number 
of stars per element of resolution. 

Our map looks ver y similar to the extinction map of 
Kai nulainen et al.l (l2006l) in terms of structure and sensitivity. 
Their map was also based on 2MASS data and a similar method 
was used to compute the extinction. The main differences are 
that they chose a higher angular resolution (2') and used only 
the 2MASS sources with a signal-to-noise ratio larger than 10 
in J, H, and K, at the cost of a few dozens of pixels being 
empty due to the lack of background stars. The extinction maps 
of Cambresv et al. (1997|) and|Cambresv (1999), based on star 
counts in the J and R bands, respectively, are more sensitive in 
the low-extinction regions than our map, but they miss the high - 
extinction regions (see also discussion in lCambresv et al.ll2002h . 

2.3. Archival C"*0 1-0 observations witti SEST 

The eastern part of Cha I was mapped with the Swedish 
ESQ Submill i meter Telescope (SEST) in the C"*0 1-0 line by 
Hai kala et al ]('2005). The angular resolution is 45" but the map 
is highly undersampled since the observations were done with 
a step of r. The median rms noise level is 0.1 K for a chan- 
nel width of 43 kHz (0.12 km s"'), but occasionally goes up 
to 0.22 K at certain positions. We retrieved the publicly avail- 
able FITS data cube from the A&A website. We reprojected 
the cube from B1950 to J2000 equatorial coordinates in radio 
projection with a projection center at (a, 5)j2oo{) - (11*'01'"24\ 
-77°15'00")- We applied a velocity correction of -1-0.25 km s"' 
to all spectra (as recommended by L. Haikala, pri v. comm.). 

3. LABOCA data reduction 

3.1. Reduction method 

The LABOCA data were reduced with the BoA softwar fol- 
lowing the procedures described in Sect. 10.2 of ISiringo et al.l 
(120091) and Sect. 3.1 of SchullergtaD (l2009l) . The removal of 
the correlated noise was done with the median noise method ap- 
plied to all pixels with 5 iterations and a relative gain of 0.9. 
Subsequently, the correlated noise computed separately for each 
group of pixels belonging to the same amplifier box was re- 
moved (BoA function correlbox with 3 iterations and a relative 
gain of 0.95), as well as the correlated noise computed sepa- 
rately for each group of pixels connected to the same read-out 
cable (BoA function correlgroup with 3 iterations and a relative 
gain of 0.8). 

The power at frequencies below 0. 1 Hz in the Fourier domain 
was partially filtered out to reduce the 1 // noise. It was replaced 
with the average power measured between 0. 1 and 0. 15 Hz (BoA 
function flattenFreq). Since the OTF scans were performed with 
a mapping speed of 2 arcmin s the frequency 0.1 Hz corre- 
sponds to an angular scale of 20'. The analysis is diff'erent for 
the spiral scans. An individual spiral subscan had a duration of 
40 s, an angular speed of 60 deg s and a radius linearly in- 
creasing from 2' to 3'. A frequency of 0.1 Hz corresponds to 
10 s, i.e. about 1.7 turns. As a result, the typical angular scale 

^ see http:// www.mpifr-bonn.mpg de/div/submmtech/software/boa/| 
boa_main.html. 
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associated with the cutoff frequency of 0.1 Hz is the spiral "di- 
ameter" which varied between 4.5' and 6'. In addition to this 
low-frequency "flattening", a first order baseline was subtracted 
scan-wise. The gridding was done with a cell size of 6.1" and the 
map was smoothed with a Gaussian kernel of size 9" (FWHM). 
The angular resolution of the final map is 21.2" (HPBW) and 
the rms noise level is 12 mJy/21.2"-beam (see Sect. 14. lb . 

This whole process was performed 21 times in an iterative 
way. The pixel values in the map produced at iteration (resp. 
1) were set to zero below a signal-to-noise ratio of 4. The loca- 
tion of the remaining signal was used as an astronomical source 
model to mask the raw data at the start of iteration 1 (resp. 2). 
The mask was removed at the end of iteration 1 (resp. 2) to com- 
pute the reduced map. Iterations 3 to 20 were performed differ- 
ently: the pixel values of the map resulting from the previous 
iteration were set to zero below a signal-to-noise ratio of 2.5. 
The remaining signal was subtracted from the raw data before 
reduction and added back after reduction. Masking or subtract- 
ing the significant astronomical signal found after one iteration 
permits to protect these regions at the next iteration when re- 
moving the correlated signal, flattening the low-frequency part 
of the power spectrum, and subtracting the first-order baseline. 
In this way, negative bowls around strong sources are much re- 
duced and more extended emission can be recovered. 




-50 50 100 

Flux density (mJy/21.2"— beam) 



Fig. 3. Flux density distribution in the full map of Cha I (up- 
per histogram), in fields Cha-Center and Cha-South (middle 
histogram), and in fields Cha-North and Cha-West (lower his- 
togram). The middle and lower histograms were shifted verti- 
cally by -500 and -1000, respectively, for clarity. A Gaussian fit 
is overlaid as a thin line on each histogram. The Gaussian stan- 
dard deviation is 12.2, 13.3, and 11.2 mJy/21.2"-beam for the 
full map, for fields Cha-Center and Cha-South, and for fields 
Cha-North and Cha-West, respectively. 



3.2. Spatial filtering and convergence 

Data reduction features dealing with spatial filtering and conver- 
gence are analysed in detail in AppendixlAl Here we give a brief 
summary. The correlated noise removal severely limits the ex- 
tent of the 870 jum emission that can be recovered (see Fig. lA.lb . 
Elongated structures are better recovered than circular ones (see 
Table lA.ll . The recovery of extended structures is improved by 
increasing the number of iterations of the data reduction. With 
21 iterations, convergence is reached for most relevant structures 
(see Figs. lA.2l and lA.3l l. 



4. Basic results 

The main assumptions made in this and the next sections to de- 
rive the physical properties of the detected sources are detailed in 
Appendix |B] These assumptions are not repeated in the follow- 
ing, except in the few cases where there could be an ambiguity. 

4. 1 . Maps of dust continuum emission in Cha I 

The final 870 fim continuum emission map of Cha I obtained 
with LABOCA is shown in Fig. |2] Pixels with a number of 
independent measurements ("coverage") smaller than 800 are 
masked. The resulting map contains 0.57 megapixels, corre- 
sponding to a total area of 1.6 deg^ (11.0 pc^). The mean and 
median coverages are 1472 and 1476 hits per pixel, respectively, 
with an integration time of 40 ms per hit. The noise distribution 
is fairly uniform and Gaussian (see Fig. [3]). The average noise 
level is 12.2 mJy/21.2"-beam. It is slightly higher for fields Cha- 
Center and Cha-South (13.3 mJy/21.2"-beam) and slightly lower 
for fields Cha-North and Cha-West (11.2 mJy/21.2"-beam). In 
the following, we will use 12 mJy/21.2"-beam as the typi- 
cal noise level. This translates into an H2 column density of 
1.1 X 10^' cm"^ for a dust mass opacity of 0.01 cm^ g"', and 
corresponds to a visual extinction Ay ~ 1.2 mag with Ry - 3.1 
(see the other assumptions in AppendixlBli. 



The dust continuum emission map of Cha I reveals several 
very compact sources and many spatially resolved sources. Very 
faint filamentary structures are also present, especially in fields 
Cha-Center and Cha-North. Figure |4] presents all the detected 
structures in more detail. The northern part (Cha-North) con- 
sists of a filamentary structure elongated along the south-north 
direction that ends with a prominent dense core dominated by a 
strong, compact source, Cha-MMS2 (source S3, Fig.|4^). Three 
additional very compact sources with weak emission are also 
detected to the east and west of the filament (S16, S20, S21). 
The northern part of field Cha-Center features one filamentary 
structure with a position angle PA ~ +60° east of north dom- 
inated by a prominent dense core that contains two compact 
sources and Cha-MMSl (S4, SIO, and CI, Fig. it). Three ad- 
ditional very compact sources are also detected outside the fila- 
ment, in the south-east, north-west, and north-east, respectively 
(S9, SI 1, S19). The region covered by the southern part of field 
Cha-Center and the northern part of field Cha-South is very 
clumpy with five very compact sources (SI, S2, S6, S7, S12) and 
many extended structures (Fig. |4j;). The southern part of field 
Cha-South contains a somewhat more isolated extended struc- 
ture (Fig.HJl). Finally, field Cha-West features relatively isolated 
sources: four compact sources (S5, S8, S17, S18), one extended 
structure and one filamentary structure (Figs.|4^-g). 

4.2. Masses traced with LABOCA and 2MASS 

The total 870 yum flux in the whole map of Cha I is about 1 15 Jy. 
Assuming a dust mass opacity k^jq - 0.01 cm^ g"', this trans- 
lates into a cloud mass of 61 M©. It corresponds to 5.9% of 
the total mass traced by CO in Cha I (1030 Mq, Mizuno et aL 
2001), 7.7% of the mass ti-acedby '^^CO (790 M^- iMizuno et al. 
199 9), and 27-32% of the mass traced by C'^O (190-230 Mq, 
Miz uno et al.lll999tlHaikala"etani200l . 

The extinction map shown in Fig.[T]traces larger scales than 
the 870 pm dust emission map. The median and mean extinc- 
tions over the 1.6 deg^ covered with LABOCA are 3.3 and 
4. 1 mag, respectively. Assuming an extinction to H2 column den- 
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Fig. 2. 870 fim continuum emission map of Cha I obtained with LABOCA at APEX. The projection type and center are the same 
as in Fig[T] The contour levels are a, 2a, 4a, So, I6a, and 32a, with a = 48 mJy/21.2"-beam, i.e. about 4 times the rms noise level. 
The flux density color scale is shown on the right. The field of view of LABOCA (10.7') and the angular resolution of the map 
(HPBW -21 .2") are shown in the lower right corner. The typical sizes above which the filtering due to the data reduction becomes 
significant in terms of peak flux density are also displayed, for weak (< 150 mJy/beam) and strong (> 150 mJy/beam) sources 
with small and large symbols, respectively, and for elliptical sources with aspect ratio 2.5 and circular sources with ellipses and 
circles, respectively (see Appendix lA. H and Col. 2 of Table FA.ll i. The pixel size is 6.1". The red boxes in the insert are labeled like 
Figs.|4^-g and show their limits overlaid on the first 870 fim contour. 



sity conversion factor of 9.4 x 10^" cm"^ mag"^ (for Ry = 3.1, 
see Appendix lB.3! ), we derive a total gas+dust mass of 950 M©. 
However, only 62% of this mass, i.e. 590 M©, is at Ay < 6 mag. 
With the appropriate conversion factor for Ay > 6 mag (see 
Appendix lB.3b . the remaining mass is reduced to 220 Mq, yield- 
ing a more accurate estimate of 810 Mq for the total mass of 
Cha I traced with the extinction. It is roughly consistent with the 
masses traced by CO and '^CO mentioned above. Since the latter 
masses may have been integrated on a somewhat different area, 
we consider the mass derived from the extinction map as the best 



estimate to compare with. Thus the mass traced with LABOCA 
represents about 7.5% of the cloud mass. Given that the median 
extinction is 3.3 mag in the extinction map, i.e. about 2.9 times 
the rms sensitivity achieved with LABOCA, the missing 92% 
were lost not only because of a lack of sensitivity but also be- 
cause of the spatial filtering due to the correlated noise removal 
(see Appendix lA.il ). Finally, we estimate the average density of 
free particles. We assume that the depth of the cloud along the 
line of sight is equal to the square root of its projected surface, 
i.e. 3.3 pc. This yields an average density of ~ 380 cm"^. For 
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Fig. 4. Detailed 870 /im continuum emission maps of Cha I extracted from the map shown in Fig. |2] The flux density greyscale 
is shown on the right of each panel and labeled in Jy/21.2"-beam. It has been optimized to reveal the faint emission with a better 
contrast. The angular resolution of the map is shown in the lower left corner of each panel (HPBW = 21.2"). The white plus 
symbols and ellipses show the positions, sizes (FWHM), and orientations of the Gaussian sources extracted with Gaussclumps or 
fitted with GAUSS_2D in the filtered map shown in Fig.|5^. The sources are labeled like in the first column of Tables|6]|8] and|9] "C" 
stands for starless core (or Class protostellar core) and "S" for YSO (Class I or more evolved). In the following, the contour levels 
are described with the parameter b - 36 mJy/21.2"-beam, i.e. about 3 times the rms noise level, a Field Cha-North. The contour 
levels are -b (in dotted blue), b, 2b, 3b, 4b, 5b, 6b, 9b, I2b, I2tb, 24b and 3Qb. The white, thin cross marks the SEST position of 
Cha-MMS2 (iReit^urth et alJll996h . 



the same depth, the density of free particles corresponding to the 
median visual extinction is estimated to be ~ 350 cm"-'. 

5. Source extraction and classification 

5.1. Multiresolution decomposition 

iMotte et^ (Il998l 12007) used a multiresolution program based 
on wavelet transforms to decompose their continuum maps into 
dififerent scales and better estimate the properties of condensa- 
tions embedded in larger-scale structures. We follow the same 
strategy with a different filter. Appendix ICldetails our method. 

This multiresolution decomposition was performed on the 
continuum map of Cha I. The total fluxes measured in the sum 
maps at scales 3 to 7 (see definition in Appendix|C]l are listed in 
Table [U as well as the corresponding masses. About half of the 



total flux is emitted by structures smaller than ~ 200" (FWHM), 
and only 11% by structures smaller than ~ 60". The sum map 
at scale 5 and its associated smoothed map are shown in Fig.|5] 
The sum of these two maps is strictly equal to the original map 
shown in Fig.|2] 

5.2. Source extraction witli Gaussclumps 

Gaussclumps (Stutzki & Gusten" l99"ol: iKramer et anil998h and 
Clumpfind (Willia ms et al.. ,1994j) are two of the most popular 
numerical codes used to extract sources from l arge-scale molec- 
ular line data cubes and dust continuum maps. Mookeriea et al.l 
(2004) found similar core mass distributions in the massive star 
forming region RCW 106 with both algorithms, even if there 
was no one-to-one correspondence between source positions and 
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ii^iEi^oo' ii''itf"oo'' ii'os^oo^ ii^)e°oo' 




iiha°oo' ii''irf"o(i" u'ofl^tM^ ii'^e°oo" 

a (J3000) 

Fig. 4. (continued) b Northern part of field Cha-Center. The contour levels are -b (in dotted blue), b , 2b, 3b, 4b, 5b, 6b, Sb, IQb, 
12b, Ub, 18/7, 22b, 26b, 30b, and 34b. The white, thin cross marks the SEST position of Cha-MMSl (iReipurth et al.lll996h . 




I a (JEOOO) 

Fig. 4. (continued) d Field Cha-South. The contour levels are -b 
(in dotted blue), b, 2b, and 3b. 



masses. O n the other hand, depe nding on the source extraction 
algorithm, ICurtis & Richej (12010 ) recently obtained significant 
differences between the core mass distributions derived for the 
Perseus molecular cloud, as well as opposite conclusions con- 
cerning the size of protostellar versus starless cores. They men- 
tioned in their conclusions that Gaussclumps may be more re- 
liable than Clumpfind to disentangle blended sources in highly- 




a (JSDDO) 



Fig. 4. (continued) e South-eastern part of field Cha-West. The 
contour levels are -b (in dotted blue), b, 2b, 4b, and lb. 



clustered regions. Here, we decided to use Gaussclumps to iden- 
tify and extract sources in the continuum map of Cha I. 
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ii''io"\H)' ii''os°'oor 11^)6^0' n'^"\H)' ii''o3:^or ll^)oroo' 




u^nToor ii^teToo* ii''os"oor ii*04"m" ii''os°oor 



Fig. 4. (continued) c Southern part of field Cha-Center. The contour levels are -b (in dotted blue), b, 2b, 3b, 4b, 5b, 6b, 9b, I5b, 
2lb, 32b, 48/7, and 64b. 



Table 1. Continuum flux distribution in Cha I. 



Scale 


Typical size 


Flux 


% of total flux 


Mass 






(Jy) 




(Mo) 


(1) 


(2) 


(3) 


(4) 


(5) 


3 


< 60" 


12.3 


11 


6.5 


4 


< 120" 


25.8 


22 


13.6 


5 


< 200" 


56.3 


49 


29.8 


6 


< 300" 


94.4 


82 


49.9 


7 


-all 


113.9 


99 


60.3 




all 


114.9 


100 


60.8 



Notes. The mass is computed assuming a dust mass opacity Kg-jo = 
0.01 cm^ g"'. The last row corresponds to the full map, while rows 
1 to 5 correspond to the sum of the filtered maps up to scale ; listed in 
the first column (i.e. the sum map at scale Column 2 gives the range 
of sizes of the sources that contribute significantly to the emission with 
more than 40% of their peak flux density (see Col. 2 of Table lC.lT l. 



We se t all three stiffness parameters to 1, as recommended 
bv lKramer et al. ( 1998). The initial guesses for the aperture cut- 
off and aperture FWHM were set to 8 and 3 times the angular 
resolution (HPBW), respectively. The initial guess for the source 



FWHM was set to 1 .SxHPBW, and the peak flux density thresh- 
old to 60 mJy/21.2"-beam, i.e. 5cr to secure the detections. 

Gaussclumps was applied to all sum maps. The sum maps 
at scales 1 to 7 were decomposed into 10, 18, 28, 42, 84, 103, 
and 114 Gaussian sources, respectively, while the full map was 
decomposed into 1 16 Gaussian sources. These counts do not in- 
clude the sources found too close to the noisier map edges (cov- 
erage < 1250), which we consider as artefacts. 

5.3. Source classification 

We now consider the results obtained with Gaussclumps for the 
sum map at scale 5 (i.e. the map shown in Fig. |5^), which is 
a good scale to characterize sources with FWHM < 120" as 
shown in Appendix ICl and Table ICTI The positions, sizes, ori- 
entations, and indexes of the 84 extracted Gaussian sources are 
listed in Table |2]in the order in which Gaussclumps found them. 
We looked for associations with sources in the SIMBAD astro- 
nomical database that provides basic data, cross-identifications, 
bibliography and measurements for astronomical objects outside 
the solar system. We used SIMBAD4 (release 1.148) as of April 
19th, 2010, and removed the objects that did not correspond to 
young stellar objects or stars. The SIMBAD sources that are lo- 
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Table 2. Sources extracted with Gaussclumps in the 870 fxm continuum sum map of Cha I at scale 5, and possible associations 
found in the SIMBAD database. 



Afgci 




R.A. 


Decl. 




/lot" 


maj." 


min.° 


P.A." 


— 777; 

S" 


Typ." 


SIMBAD possible 


Dist." 






(J200()) 


(J2()00) 


(Jy/beam) 


(Jy) 


(") 


(") 


(°) 


(.") 




associations'' 


(") 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


11 


:08:03.22 


-77:39:17.9 


2.580 


3.093 


23.5 


22.9 


50.0 


23.2 


S 


V* CU Cha (Or*) 


0.5 
























BRAN 34 IB (BNe) 


1.3 


2 


11 


:08:38.96 


-77:43:52.9 


1.206 


1.893 


26.6 


26.6 


-33.2 


26.6 


s 


Cha IR NEBULA (TT*) 


6.2 


3 


11 


:10:00.13 


-76:34:59.1 


1.169 


1.501 


25.1 


23.0 


61.3 


24.0 


s 


V* WW Cha (Or*) 


1.2 
























CHXR 44 (X) 


6.6 
























[RNC96] Cha-MMS 2 (mm) 


7.5 


4 


11 


:06:31.94 


-77:23:38.9 


1.053 


6.350 


55.1 


49.3 


41.5 


52.1 


c 


[RNC96] Cha-MMS la (mm) 


6.9 


5 


11 


:06:46.44 


-77:22:32.2 


0.607 


2.418 


57.8 


31.0 


80.9 


42.3 


s 


[PMG2001] NIR 69 (IR) 


2.2 
























IRAS 11051-7706 (Y*0) 


2.9 
























[PMG2001] NIR 72 (IR) 


5.9 
























[PMG2001] NIR 74 (IR) 


6.4 
























[PMG2001] NIR 80 (IR) 


12.8 


6 


11 


:06:15.51 


-77:24:04.9 


0.306 


3.067 


98.2 


45.9 


-7.0 


67.2 


c 


[PMG2001] NIR 16 (IR) 


38.2 


7 


10 


:58: 16.92 


-77:17:18.3 


0.293 


0.314 


22.7 


21.3 


-34.0 


21.9 


s 


V* SZ Cha (Or*) 


1.3 


8 


11 


:07:01.85 


-77:23:00.8 


0.283 


2.197 


76.1 


45.9 


53.7 


59.1 


c 


ISO-Chal 86 (Y*0) 


14.6 
























IRAS 11057-7706 (Y*0) 


26.4 
























[PCW91] Ced 110 IRS 6A (IR) 


27.0 


9 


11 


: 04: 16.41 


-77:47:32.1 


0.209 


1.249 


53.4 


50.4 


-78.1 


51.9 


c 






10 


11 


:02:24.98 


-77:33:36.3 


0.195 


0.197 


21.4 


21.2 


-63.0 


21.3 


s 


V* CS Cha (Or*) 


0.6 


11 


11 


:08:50.40 


-77:43:46.7 


0.167 


0.773 


60.3 


34.5 


5.2 


45.6 


c 


HJMC 9-1 (IR) 


8.2 


12 


11 


:08: 15.28 


-77:33:51.9 


0.168 


0.173 


21.9 


21.2 


-63.0 


21.6 


s 


Glass la (TT*) 


1.2 
























SSTc2d 11 10815.3-773353 (*) 


1.4 


13 


11 


:08:00.07 


-77:38:43.8 


0.161 


0.412 


44.3 


26.0 


79.7 


33.9 


c 


HH 920B (HH) 


6.8 
























2MJ1 1075792-7738449 (TT*) 


7.0 
























2MJ1 1080297-7738425 (Y*0) 


9.4 
























[HTN2003] HD 97048 2 (IR) 


11.5 
























ESO-HA 562 (Y*0) 


12.3 


14 


11 


:09:45.85 


-76:34:49.7 


0.148 


1.552 


98.3 


48.1 


-7.9 


68.8 


c 


2M 111094621-7634463 (TT*) 


3.5 
























Hn low (Em*) 


4.5 
























[ACA2003] Cha 14 (Em*) 


7.3 
























OTS 20 (*) 


13.0 
























[ACA2003]ChaI3(pr*) 


16.2 


15 


11 


:06:59.93 


-77:22:00.6 


0.154 


0.276 


31.0 


26.1 


15.0 


28.4 


c 


[PMG2001] NIR 103 (IR) 


7.5 
























[PMG2001] NIR 96 (IR) 


7.9 
























[PMG2001] NIR 101 (IR) 


10.8 
























[PMG2001] NIR 93 (IR) 


11.9 


16 


10 


:56:42.95 


-77:04:05.1 


0.134 


1.278 


82.6 


51.8 


23.0 


65.4 


c 






17 


10 


:56:30.56 


-77:11:39.8 


0.137 


0.142 


22.0 


21.2 


-88.6 


21.6 


s 


V* SY Cha (Or*) 


0.6 


18 


11 


:09:47.42 


-77:26:31.5 


0.129 


0.132 


21.8 


21.2 


-10.2 


21.5 


s 


2M 11 1094742-7726290 (Y*0) 


2.4 


19 


11 


:06:39.88 


-77:24:18.1 


0.130 


0.521 


68.1 


26.5 


53.7 


42.4 


R 






20 


11 


:07:21.34 


-77:22:12.1 


0.119 


0.121 


21.6 


21.2 


-5.7 


21.4 


s 


DENIS-P Jl 107.4-7722 (Y*0) 


1.3 
























[PMG2001] NIR 137 (IR) 


3.8 
























[PMG2001] NIR 140 (IR) 


4.4 


21 


11 


:06:06.48 


-77:25:06.7 


0.116 


0.583 


65.9 


34.3 


-20.9 


47.5 


c 






22 


11 


:04:23.35 


-77:18:07.8 


0.114 


0.115 


21.4 


21.2 


-19.0 


21.3 


s 


HH 48 IRS (HH) 


2.0 
























HH 48E (HH) 


7.5 


23 


11 


:02: 10.85 


-77:42:31.8 


0.111 


1.434 


126.7 


45.8 


38.6 


76.2 


c 






24 


11 


:09:26.24 


-76:33:28.3 


0.111 


0.455 


56.3 


32.7 


27.1 


42.9 


c 


DENIS-P Jl 109.4-7633 (Y*0) 


5.2 
























[ACA2003]ChaI2(pr*) 


8.1 


25 


11 


:03:36.59 


-77:59:19.3 


0.104 


1.023 


69.4 


63.7 


85.8 


66.5 


c 






26 


11 


:07:17.13 


-77:23:20.2 


0.109 


0.403 


49.8 


33.6 


-17.3 


40.9 


c 


[PMG2001] NIR 134 (IR) 


4.8 
























DENIS-P Jl 107.3-7723 (Y*0) 


11.9 
























[PMG2001] NIR 135 (IR) 


15.4 
























[PMG2001] NIR 142 (IR) 


16.7 


27 


11 


:06:35.35 


-77:22:52.8 


0.106 


0.263 


52.3 


21.4 


-74.9 


33.5 


R 


lAXGJl 10700-7728 (X) 


25.3 


28 


11 


:07:43.89 


-77:39:41.6 


0.093 


0.094 


21.4 


21.2 


-35.9 


21.3 


s 


Ass Cha T 2-28 (TT*) 


0.8 


29 


11 


:10:01.05 


-76:36:32.9 


0.090 


0.999 


102.4 


48.6 


-8.0 


70.5 


c 


OTS 37 (*) 


28.8 


30 


11 


:06:21.67 


-77:43:58.4 


0.089 


0.392 


47.4 


41.7 


-60.9 


44.5 


c 






31 


11 


: 10: 12.26 


-76:35:30.1 


0.093 


0.211 


35.8 


28.5 


-79.5 


32.0 


Sc 


2M Jl 1101 141-7635292 (pr*) 


3.0 


32 


11 


:03:21.97 


-77:36:04.4 


0.088 


0.554 


67.3 


42.4 


-77.0 


53.4 


c 




33 


11 


:04:49.58 


-77:45:53.5 


0.086 


0.926 


102.1 


47.7 


75.1 


69.8 


c 
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Table 2. continued. 







R.A. 


Decl. 


/peak 


ftot^ 


maj." 


min." 


P.A." 




Typ.'- 


SIMBAD possible 


Dist.'' 






(J2000) 


02000) 


(Jy/beam) 


(Jy) 


(") 


(") 


(°) 


(") 




associations 


(") 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


34 


11 


03:39.07 


-77:47:54.0 


0.083 


0.529 


76.3 


37.5 


60.7 


53.5 


C 






35 


11 


11:27.78 


-77:15:45.2 


0.078 


0.507 


73.5 


39.9 


83.3 


54.1 


C 






36 


11 


04:07.89 


-77:48:13.1 


0.082 


0.197 


37.4 


29.0 


-12.9 


32.9 


c 






37 


11 


04:01.79 


-77:58:50.4 


0.078 


1.016 


90.7 


64.3 


-8.9 


76.3 


c 






38 


11 


03:04.43 


-77:34:19.4 


0.079 


0.262 


40.5 


37.0 


-62.7 


38.7 


c 






39 


11 


07:17.30 


-77:38:42.3 


0.077 


0.307 


50.3 


35.8 


-10.8 


42.5 


c 


HH 927 (HH) 


12.3 


40 


11 


09:03.86 


-77:43:34.6 


0.081 


0.346 


51.4 


37.2 


-9.1 


43.8 


c 






41 


11 


02:10.46 


-77:37:22.9 


0.077 


0.830 


111.5 


43.7 


-10.3 


69.8 


c 






42 


11 


09:39.48 


-76:36:04.6 


0.080 


0.400 


61.1 


36.7 


-83.9 


47.3 


c 






43 


11 


07:47.07 


-77:35:46.8 


0.075 


0.319 


59.6 


32.1 


22.0 


43.7 


c 






44 


11 


05:57.64 


-77:24:18.9 


0.082 


0.323 


49.3 


36.1 


44.8 


42.2 


c 






45 


11 


00:43.98 


-77:30:38.5 


0.072 


0.862 


105.2 


51.2 


-80.5 


73.4 


c 






46 


11 


00:58.29 


-77:29:41.2 


0.074 


0.528 


71.6 


44.6 


-33.4 


56.5 


c 






47 


11 


07:32.72 


-77:36:23.0 


0.075 


0.752 


91.5 


49.5 


-61.2 


67.3 


c 






48 


11 


10:01.60 


-76:34:54.3 


0.077 


0.078 


21.4 


21.2 


-6.3 


21.3 


R 


CHXR 44 (X) 
V* WW Cha (Or*) 


1.8 
6.3 


49 


11 


02:33.40 


-77:42:20.1 


0.075 


0.348 


58.2 


36.0 


0.0 


45.8 


c 






50 


10 


50:56.54 


-77:04:29.1 


0.072 


0.952 


139.3 


42.7 


-75.1 


77.1 


c 






51 


11 


08:12.67 


-77:39:12.1 


0.080 


0.130 


34.6 


21.2 


0.3 


27.1 


R 






52 


11 


09:15.17 


-77:18:10.7 


0.071 


0.235 


55.8 


26.6 


38.4 


38.5 


C 






53 


11 


06:23.17 


-77:41:42.7 


0.074 


0.271 


49.1 


33.7 


79.1 


40.7 


c 






54 


11 


09:56.09 


-77:14:16.7 


0.070 


0.365 


77.6 


30.1 


86.9 


48.3 


c 






55 


11 


02:10.12 


-77:39:29.2 


0.071 


0.587 


83.4 


44.6 


36.2 


61.0 


c 






56 


11 


10:03.64 


-76:38:02.6 


0.068 


0.445 


82.3 


35.8 


-46.8 


54.2 


c 






57 


11 


06:03.63 


-77:34:52.7 


0.070 


0.426 


58.7 


46.6 


-81.7 


52.3 


c 






58 


11 


02:58.44 


-77:38:42.5 


0.070 


0.241 


48.3 


31.9 


-0.7 


39.3 


c 






59 


11 


07:15.99 


-77:24:35.6 


0.068 


0.192 


48.7 


26.0 


44.0 


35.6 


c 


[PMG2001] NIR 129 (IR) 
[PMG2001] NIR 130 (IR) 
[PMG2001] NIR 126 (IR) 
[PMG2001] NIR 136 (IR) 


8.7 
8.8 
12.8 
14.0 


60 


11 


09:51.10 


-76:51:19.2 


0.066 


0.598 


90.0 


45.6 


34.2 


64.0 


c 


ISO-Chal 206 (IR) 


11.7 


61 


11 


08:59.13 


-77:45:08.6 


0.067 


0.292 


46.5 


42.4 


-76.4 


44.4 


c 






62 


11 


09:53.92 


-76:34:23.2 


0.070 


0.071 


21.5 


21.2 


10.2 


21.3 


Sc 


HH 914 (HH) 

CED 112 IRS 4 (TT*) 


1.8 
2.9 


63 


11 


10:01.87 


-76:33:11.5 


0.068 


0.199 


54.5 


24.3 


-33.5 


36.4 


Sc 


2M Jl 1100336-76331 11 (Y*0) 
[FL2004] 55 (pr*) 


5.3 
18.8 


64 


11 


11:22.58 


-77:16:24.4 


0.066 


0.229 


73.2 


21.4 


-87.1 


39.6 


c 






65 


11 


09:00.36 


-77:40:12.6 


0.065 


0.112 


36.6 


21.2 


-66.1 


27.9 


c 






66 


11 


05:38.45 


-77:49:12.5 


0.065 


0.069 


22.4 


21.4 


69.9 


21.9 


A 






67 


11 


07:24.50 


-77:39:40.6 


0.065 


0.459 


87.6 


36.2 


42.5 


56.3 


c 






68 


11 


02:33.46 


-78:01:38.2 


0.063 


0.170 


50.7 


24.0 


87.8 


34.9 


c 






69 


11 


06:24.83 


-77:22:58.5 


0.070 


0.112 


33.8 


21.2 


53.0 


26.8 


R 


[PMG2001] NIR 35 (IR) 


5.3 


70 


11 


04:47.42 


-77:44:58.1 


0.064 


0.154 


40.2 


26.9 


-63.6 


32.9 


C 






71 


11 


09:07.16 


-77:23:55.4 


0.063 


0.457 


60.6 


54.3 


52.1 


57.4 


C 






72 


11 


07:35.16 


-77:21:45.8 


0.064 


0.158 


34.9 


31.9 


-86.8 


33.4 


C 






73 


11 


08:03.98 


-77:20:53.1 


0.061 


0.188 


50.5 


27.8 


-55.9 


37.4 


C 






74 


11 


10:23.04 


-76:38:40.4 


0.063 


0.295 


52.8 


40.1 


13.1 


46.0 


C 






75 


11 


06:30.97 


-77:24:25.2 


0.068 


0.090 


27.7 


21.4 


66.2 


24.4 


R 






76 


11 


08:45.48 


-77:44:22.7 


0.063 


0.080 


27.1 


21.2 


68.3 


24.0 


R 






"7"? 


1 1 




- /o:4/:UU.2 


U.UOZ 


U.U/Z 


24.6 


TIT 

zi.z 


on c 

zy.3 


ll>j 


c 
J) 


2M J i 1 i 1 iUoi- /o4i3 /4 (I'^O) 


J. 3 


78 


11 


07:28.55 


-77:19:40.6 


0.060 


0.145 


36.6 


29.9 


-70.1 


33.1 


c 






79 


11 


06:18.59 


-77:44:40.1 


0.061 


0.159 


40.7 


28.8 


87.1 


34.2 


c 






80 


11 


05:28.19 


-77:41:18.1 


0.060 


0.137 


46.3 


22.4 


-23.4 


32.2 


c 






81 


10 


56:49.93 


-77:03:02.8 


0.061 


0.280 


51.1 


40.1 


-86.3 


45.3 


c 






82 


11 


02:12.48 


-77:41:19.1 


0.063 


0.214 


48.5 


31.6 


82.2 


39.2 


c 






83 


11 


03:00.75 


-77:49:18.2 


0.061 


0.118 


41.0 


21.4 


-82.5 


29.6 


c 






84 


11 


09:53.67 


-76:42:12.0 


0.058 


0.350 


61.6 


43.7 


23.5 


51.9 


c 







Notes. Peak flux density (in Jy/21.2"-beam), total flux, FWHM along the major and minor axes, and position angle (east from north) of the 
fitted Gaussian. Mean source size, equal to the geometrical mean of the major and minor FWHM. Type of source based on the associations 
found in the SIMBAD database (see Tabled. Source found in the SIMBAD database within the FWHM ellipse, its type in parentheses, and its 
distance to the peak position. 2MASS was abbreviated to 2M. The abbreviated references in square brackets and the source types in parentheses 
are defined in SIMBAD. 
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ii^iiToo* ii^wPtoo* ic^oHw' ii^iiToo* ic^oHw' 




ii^nnjoT iiNMrooT loi'iMrtH" ii^nnjor uNMnmr loi'iMrtH" 



Fig. 5. a 870 /vm continuum emission sum map of Cha I at scale 5 (see Sect. 15.11 ). The contour levels are -a (in dashed blue), a, 2a, 
4a, 8fl, I6a, and 32a, with a = 48 mJy/21.2"-beam, i.e. about 4 times the rms noise level, b Smoothed map, i.e. residuals, at scale 
5. The contour levels are -c (in dashed blue), c, 2c, 4c, and 8c, with c = 13.5 mJy/21.2"-beam, i.e. about 4.5 times the rms noise 
level in this map. The labels indicate the filaments listed in Table [TOl The greyscales of both maps are different. The sum of these 
two maps is strictly equal to the original map (Fig.|2]l. 




Fig. 4. (continued) f North-eastern part of field Cha- West. The 
contour levels are -b (in dotted blue), b, 2b, 3b, and 4b. 



cated within the FWHM ellipse of each Gaussclumps source 
after this filtering are listed in Table |2] along with their type and 
their distance to the fitted peak position. Based on these possible 




Fig. 4. (continued) g North-western part of field Cha-West. The 
contour levels are -b (in dotted blue), b, and 2b. 

associations, we classify the sources found with Gaussclumps 
into five categories (Col. 10 of Table|2]i: 

- S: very compact or unresolved source associated with a 
young stellar object (Class I or more evolved) the position 
of which agrees within a few arcsec (13 sources), 

- Sc: tentative association like S, but with some uncertainty 
due to the somewhat extended 870 fim emission (3 sources), 

- R: source that looks like a residual (departure from 
Gaussianity) of a stronger nearby source and may therefore 
not be real (7 sources), 

- A: likely artefact, i.e. unresolved source with no SIMBAD 
association (1 source). 
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Table 3. Classification of the sources extracted with 
Gaussclumps from the 870 /vm continuum sum map of Cha 1 
at scale 5. 



Category 


Number 


Definition 


(1) 


(2) 


(3) 


C 


60 


starless core or Class protostar 


s 


13 


young stellar object (Class 1 or more evolved) 


Sc 


3 


candidate association with a YSO 


R 


7 


residual of a stronger nearby source 


A 


1 


artefact 



- C: source that looks like a bona-fide starless core or Class 
protostar (60 sources). 

This classification is summarized in Table[3] 

Twelv^ LABOCA sources classified as bona-fide starless 
cores (C) in Table|2]have a possible SIMBAD association within 
their FWHM ellipse {N^d = 6, 8, 11, 13 14, 15, 24, 26, 29, 39, 
59, and 60). They are still classified as starless cores because ei- 
ther they are extended and do not show any evidence for a com- 
pact structure at the position of the possibly associated SIMBAD 
source, or their peak position is significantly offset from it. We 
consider that the SIMBAD source is either not physically as- 
sociated with the dense core (chance association, especially in 
the large dense core in Cha-Center where an IR cluster overlaps 
with its outskirts) or it is embedded in the dense core but the 
latter may still form an additional star However, we cannot ex- 
clude that a fraction of these twelve dense cores are remnant of 
a previous episode of star formation and will never form new 
stars. 

Cha I was observed with the Spitzer instruments IRAC 
and MI PS in the fram ework of the Gould's Belt legacy pro- 
gram (.Allen et al.1 12006) . The maps are currently being anal- 
ysed (J0rgensen et al., in prep.) and were used to search for 
possible associations with the LABOCA sources. Among the 
60 LABOCA sources classified as bona-fide starless cores (C), 
7 have a compact counterpart in emission in the 24 or 70 yum 
Spitzer MIPS maps within a radius of 15". These Spitzer sources 
are listed in Table H) The first one (A^sp = 1) is associated with 
source A^gd - 4 (C ha-MMSl) and was already reported by 
[Belloche"etal]( l2006l) . who argued that it is a very young Class 
protostar, or maybe even at the stage of the first hydrostatic core. 
At this stage, the mass in the protostellar envelope completely 
dominates the mass of the stellar embryo. Therefore it makes 
sense to compare its properties with those of starless cores and 
we keep it in category C. Note that the peak of the 870 jum emis- 
sion is actually a bit shifted compared to the position fitted by 
Gaussclumps and reported in Table|2] The deviation comes from 
the asymmetry of the core. By eye, we measure a peak posi- 
tion {a,5)n(.m = (ll''06'"32^.94, -77°23'35'.'2), i.e. much closer 
to the Spitzer position with a relative offset (-0.6", -1.1"), which 
strengthens the association of the Spitzer source with the 870 
source A^gd - 4 even further 

The second and third Spitzer sources (A^sp - 2 and 3) are 
coincident within 0.3" with SIMBAD YSOs already listed in 
Table |2l Their mid-IR slope indicates that they are Class II ob- 
jects (J0rgensen et al., in prep.). The fourth and fifth Spitzer 
sources (A^sp - 4 and 5) are also coincident with SIMBAD ob- 
jects, but only within 3", which makes the association less se- 
cure. Their mid-IR slope indicates that they are Class II and 

* We do not count Cha-MMSla since it is only a 1.3 mm source in 
SIMBAD. 



flat-spectrum objects, respectively. The four 870 ;um sources 
A^gd = 13, 14, 24, and 26 possibly associated with the Spitzer 
sources A^sp = 2 to 5 are extended. Since they do not show ev- 
idence for a compact structure toward the Spitzer source, we 
still classify them as bona-fide starless cores (C). The last two 
Spitzer sources (A^sp - 6 and 7) are Class I and II objects, re- 
spectively, according to their mid-IR slope, but they do not have 
any SIMBAD counterpart. They are located ~ 6" and ~ 7" from 
the peak of the 870 jjm sources A^gd = 37 and 39, respectively. 
Since the latter are also extended and do not show evidence for 
a compact structure, we also still classify them as bona-fide star- 
less cores (C). 

In summary, out of 84 sources found with Gaussclumps in 
the 870 //m map of Cha I, 76 look like real sources. Sixteen of 
these sources (21%) are associated with a young stellar object 
(Class I or more evolved) and most likely trace a circumstel- 
lar disk and/or a residual circumstellar envelope. The remaining 
60 sources (79%) appear to be starless based on Spitzer or con- 
tain an embedded Class protostar or first hydrostatic core (one 
case). The position, sizes, and orientation of these 76 sources, 
plus the additional sources mentioned in Sect. 15.41 below, are 
shown in Figs. |4^ to g. The labels correspond to the first col- 
umn of Tables |6]|8] and|9] "C" stands for starless core (or Class 
protostellar core) and "S" for YSO (Class I or more evolved). 

5.4. Additional sources 

Since the detection threshold for Gaussclumps was set to 5cr, 
compact sources with a SIMBAD counterpart and a 870 fim de- 
tection between 3.5cr and ~ 5cr that may be trusted based on 
the association were missed by Gaussclumps. Therefore, we also 
looked for 870 pm emission above 3.5cr in the sum map at scale 
5 at the position of each source in the SIMBAD database. Clear 
associations between a SIMBAD source and a compact 870 pm 
emission were visually selected and all SIMBAD sources as- 
sociated with extended emission but no clear peak at 870 pm 
were discarded, in particular those in the prominent dense cores 
around Cha-MMSl and Cha-MMS2. As a result, there are five 
additional compact 870 pm sources with a likely SIMBAD as- 
sociation and a formal peak flux density above 3.5cr. They are 
listed in Table |5] along with the parameters derived from a 
Gaussian fit performed with the task GAUSS_2D in GRECQ 
The fitted peak flux density is below 3.5cr in one case. Given 
the number of SIMBAD sources in the Cha I field (1045 ob- 
jects over 1.6 deg^), the probability of chance association within 
a radius of 5" is 0.4% only, so we are confident that these five 
additional SIMBAD associations are real. 



6. Analysis 

6. 1 . Starless cores 

The properties of the 60 starless (or Class 0) sources are listed 
in Tableland their distribution is shown in Figs. |6] and [8] The 
column density (Col. 6) and masses (Cols. 8, 9, and 10) are com- 
puted with the fluxes fitted with Gaussclumps or directly mea- 
sured in the sum map at scale 5, assuming a dust mass opacity 
'<'870 - 0.01 cm^ g"'. As a caveat, we remind that the assumption 
of a uniform temperature may be inaccurate and bias the mea- 
surements of the masses and column densities, as well as the 
mass concentration (or equivalently the density contrast). Since 

' GREG is part of the collection of softwares GILDAS (see 
|http://www.iram .fr/lRAM FR/GlLDAS^ . 
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Table 4. Properties of the Spitzer sources possibly associated with a 870 jim source classified C in Table |2] 



A'sp 


R.A. 


Decl. 


a" 




Dist.^ 


SIMBAD possible 


Dist."* 




(J2000) 


(J2000) 






(") 


associations'' 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


11:06:33.13 


-77:23:35.1 


- 


4 


5.5 






2 


11:07:57.96 


-77:38:45.2 


-0.3 


13 


7.0 


2MJ1 1075792-7738449 (TT*) 


0.3 


3 


11:09:46.19 


-76:34:46.6 


-0.4 


14 


3.3 


2M J 1109462 1-7634463 (TT*) 


0.2 














Hn low (Em*) 


1.6 














[ACA2003] Clia 14 (Em*) 


4.8 


4 


11:09:26.09 


-76:33:33.9 


-1.0 


24 


5.7 


DENIS-PJI 109.4-7633 (Y*0) 


3.0 


5 


11:07:16.19 


-77:23:06.9 


-0.1 


26 


13.5 


DENIS-PJl 107.3-7723 (Y*0) 


3.0 


6 


11:04:01.15 


-77:58:44.7 


1.0 


37 


6.1 






7 


11:07:16.32 


-77:38:48.3 


-0.3 


39 


6.9 







Notes. Mid-IR spectral slope (j0rgensen et al., in prep.). A dash means that it is not detected shortward of 24 ^m. Possibly associated 
Gaussclumps 870 pm source (same number as in Table |2]l. Distance to the peak of the Gaussclumps 870 pm source. Source found in the 
SIMBAD database, its type in parentheses, and its distance to the Spitzer source. Only sources with a distance smaller than 5" are listed. 2MASS 
was abbreviated to 2M. 

Table 5. Additional compact sources with SIMBAD association with a formal peak flux density above 3.5cr in the 870 //m contin- 
uum map of Cha I filtered up to scale 5. 



A'add 


R.A. 


Decl. 


J peak 


/lot" 


maj." 


min." 


PA." 


5* 


SIMBAD possible 


Dist.'^ 




(J2000) 


(J2000) 


(Jy/beam) 


(Jy) 


(") 


(") 


(°) 


(") 


associations'^ 


(") 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


10:59:01.28 


-77:22:40.4 


0.055 


0.062 


24.4 


20.7 


17.9 


22.5 


V* TW Cha (Or*) 


0.7 


2 


10:55:59.55 


-77:24:35.9 


0.032 


0.075 


56.9 


18.6 


69.5 


32.5 


Ass Cha T 2-3 B (Y*0) 


4.1 




















2M J10555973-7724399 (Y*0) 


4.1 




















V* SX Cha (Or*) 


4.2 


3 


11:10:49.76 


-77:17:51.0 


0.047 


0.045 


23.2 


18.4 


-22.5 


20.7 


HBC 584 (Em*) 


0.9 


4 


11:06:25.43 


-76:33:44.8 


0.059 


0.044 


20.4 


16.6 


-88.3 


18.4 


2MJ1 1062554-7633418 (Y*0) 


3.0 


5 


11:12:28.26 


-76:44:24.2 


0.061 


0.041 


18.6 


16.3 


87.2 


17.4 


V* CV Cha (Or*) 


2.7 




















BRAN341D (RNe) 


4.0 




















2E 1110.9-7628 (X) 


5.5 



Notes. Peak flux density (in Jy/21.2"-beam), total flux, FWHM along the major and minor axes, and position angle (east from north) of the 
fitted Gaussian. The fitted peak flux density may slightly depart from the formal peak flux density measured directly in the map. Mean source 
size, equal to the geometrical mean of the major and minor FWHM. Source found in the SIMBAD database within the FWHM ellipse, its type 
in parentheses, and its distance to the peak position. 2MASS was abbreviated to 2M. 



only one of the 60 sources has an embedded YSO (as traced at 
24 fim with Spitzer, see Sect. I5.3l l. the other 59 have no central 
heating and their temperature should be rather uniform. However 
a dust temperature drop toward the center of starless dense cores 
is possible (see Appendix[ 



6.1.1. Extinction 

The visual extinctions listed in Table |6] are extracted from the 
extinction map derived from 2MASS (see Sect. I2.2l i. Given the 
lower resolution of this map (HPBW = 3') compared to the 
870 fim map, it provides an estimate of the extinction of the en- 
vironment in which the 870 ;um sources are embedded. 

The 870 jum sources are found down to a visual extinction 
Ay ~ 5 mag (as traced with 2MASS at low ang ular resolution), 
but not below. A similar result was obtained by 'Johnstone et al.' 
(2004) in Ophiuchus (threshold at Ay - 7 mag), by Enoch et al. 
(l2006l) and iKirk e t al. (200d) in Per seus (threshold at Ay = 
5 mag), and by Andre et a P ('2010a) in Aquila (threshold at 
Ay = D niag, see Fig. 5 oflAndre et al. 2010b, hut see below for 
the peak of the distribution). In Taurus, Golds rnith et al.l (l2008l) 
also found a threshold at Ay = 6 mag based on the d istribution of 
YSOs , somewhat lower than the earlier findings of lOnishi et al.l 
(Il998h who proposed a threshold at A^h, = 8 x 10^' cm"^ 



(Ay ~ 9 mag) based on C'^O cores and IRAS sources. In both 
latter cases, however, this threshold is not sharp but rather indi- 
cates a significant increase in the probability for star formation 
to occur 

The extinction toward the only confirmed young protostar 
in Cha I (Cha-MMSl or Cha 1 -CI in Table |6]l is high, at = 
14 mag. The distribution of extinctions of the environment in 
which the starless cores are embedded has a median of 9.1 mag 
and presents a sharp peak at Ay ~ 8-9 mag (see Fig.|6^). This is 
very similar in shape to the distribution of extinctions found for 
A quila with He rschel, which peaks at Ay ~ 7.5 mag (see Fig. 5 
of lAndre et ani20T0bl) . 

Star formation is more likely at high extinction, as discussed 
above, and in Cha I the areas of high extinction are small. We 
can further quantify the relationship between the extinction and 
the location of the dense cores using the distribution of extinc- 
tions across the map (Fig.|7^). We calculate the average surface 
density of cores Ecoies in each extinction bin by dividing the num- 
ber of starless (+ Class 0) cores in that bin by the total area of 
the map (in pc^) in that extinction range. The resulting average 
surface density of cores is plotted as a function of gas column 
density Egas (as measured by A,,) as a blue histogram in Fig.|7]3. 

If we assume that each core will ultimately form a star, 
then we could adopt the extragalactic concept of the Kennicutt- 
Schmidt law and convert the density of cores to a star formation 



13 



A. Belloche et al.: The end of star formation in Cliamaeleon I ? 

Table 6. Characteristics of starless (or Class 0) sources extracted with Gaussclumps in the 870 /^m continuum map of Cha I filtered 
up to scale 5. 



Name 




FWHM" 


FWHM" 




^peak 








^50"' 


Cm' 




50 


c * 






("X") 


(1000 AU)" 




(10-' cm" 


'-) (mag) 


(Mo) 


(Mo) 


(Mo) 


(%) 


(10' cm" 


') (10* cm" 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Cha 1 -CI 


4 


55.1 X 49.3 


7.6 X 6.7 


1.1 


92 


14 


0.56 


3.36 


1.44 


39( 1) 


49.5 


97.7 


5.1(0.1) 


Chal-C2 


6 


98.2 X 45.9 


14.4 X 6.1 


2.4 


27 


9 


0.16 


1.62 


0.58 


28( 1) 


14.4 


39.4 


3.7(0.2) 


Chal-C3 


8 


76.1 X 45.9 


11.0 X 6.1 


1.8 


25 


13 


0.15 


1.16 


0.45 


33( 2) 


13.3 


30.7 


4.3(0.2) 


Chal-C4 


9 


53.4 X 50.4 


7.3 X 6.9 


1.1 


18 


17 


0.11 


0.66 


0.29 


38( 3) 


9.8 


19.9 


4.9(0.3) 


Chal-C5 


11 


60.3 X 34.5 


8.5 X 4.1 


2.1 


15 


9 


0.088 


0.41 


0.32 


28( 2) 


7.9 


21.4 


3.7(0.3) 


Chal-C6 


13 


44.3 X 26.0 


5.8 X 2.2 


2.6 


14 


8 


0.085 


0.22 


0.33 


26( 2) 


7.6 


22.3 


3.4(0.3) 


Chal-C7 


14 


98.3 X 48.1 


14.4 X 6.5 


2.2 


13 


11 


0.078 


0.82 


0.23 


34( 3) 


6.9 


15.4 


4.5(0.4) 


Chal-C8 


15 


31.0 X 26.1 


3.4 X 2.3 


1.5 


13 


14 


0.081 


0.15 


0.22 


36( 3) 


7.2 


15.2 


4.8(0.4) 


Chal-C9 


16 


82.6x51.8 


12.0 X 7.1 


1.7 


12 


7 


0.071 


0.68 


0.20 


36( 4) 


6.3 


13.4 


4.7(0.5) 


Chal-ClO 


21 


65.9 X 34.3 


9.4 X 4.0 


2.3 


10 


9 


0.061 


0.31 


0.19 


33( 4) 


5.5 


12.8 


4.3(0.5) 


Chal-Cll 


23 


126.7 X 45.8 


18.7 X 6.1 


3.1 


9.7 


9 


0.059 


0.76 


0.16 


36( 5) 


5.2 


11.1 


4.7(0.6) 


Chal-C12 


24 


56.3 X 32.7 


7.8 X 3.7 


2.1 


9.7 


10 


0.059 


0.24 


0.14 


42( 5) 


5.2 


9.6 


5.5(0.7) 


Chal-C13 


25 


69.4 X 63.7 


9.9 X 9.0 


1.1 


9.1 


9 


0.055 


0.54 


0.16 


35( 5) 


4.9 


10.5 


4.7(0.6) 


Chal-C14 


26 


49.8 X 33.6 


6.8x3.9 


1.7 


9.5 


13 


0.057 


0.21 


0.17 


34( 4) 


5.1 


11.6 


4.4(0.6) 


Chal-C15 


29 


102.4 X 48.6 


15.0 X 6.6 


2.3 


7.9 


8 


0.048 


0.53 


0.13 


37( 6) 


4.2 


8.8 


4.8(0.8) 


Chal-C16 


30 


47.4x41.7 


6.4 X 5.4 


1.2 


7.8 


10 


0.047 


0.21 


0.11 


44( 7) 


4.2 


7.3 


5.8(1.0) 


Chal-C17 


32 


67.3 X 42.4 


9.6 X 5.5 


1.7 


7.7 


11 


0.046 


0.29 


0.11 


41(7) 


4.1 


7.7 


5.3(0.9) 


Chal-C18 


33 


102.1 X 47.7 


15.0 X 6.4 


2.3 


7.5 


13 


0.045 


0.49 


0.12 


38( 6) 


4.0 


8.0 


5.0(0.8) 


Chal-C19 


34 


76.3 X 37.5 


11.0x4.6 


2.4 


7.3 


14 


0.044 


0.28 


0.098 


45( 8) 


3.9 


6.7 


5.9(1.1) 


Chal-C20 


35 


73.5 X 39.9 


10.6 X 5.1 


2.1 


6.8 


5.2 


0.041 


0.27 


0.095 


43( 8) 


3.7 


6.5 


5.7(1.1) 


Chal-C21 


36 


37.4 X 29.0 


4.6 X 3.0 


1.6 


7.1 


19 


0.043 


0.10 


0.11 


38( 7) 


3.8 


7.7 


5.0(0.9) 


Chal-C22 


37 


90.7 X 64.3 


13.2 X 9.1 


1.5 


6.9 


9 


0.041 


0.54 


0.12 


35( 6) 


3.7 


8.1 


4.5(0.8) 


Chal-C23 


38 


40.5 X 37.0 


5.2x4.5 


1.1 


6.9 


10 


0.042 


0.14 


0.087 


48( 9) 


3.7 


5.9 


6.3(1.2) 


Chal-C24 


39 


50.3 X 35.8 


6.8 X 4.3 


1.6 


6.7 


14 


0.041 


0.16 


0.081 


50(10) 


3.6 


5.5 


6.6(1.3) 


Chal-C25 


40 


51.4 X 37.2 


7.0 X 4.6 


1.5 


7.1 


8 


0.043 


0.18 


0.12 


35( 6) 


3.8 


8.3 


4.6(0.8) 


Chal-C26 


41 


111.5 X 43.7 


16.4 X 5.7 


2.9 


6.7 


7 


0.041 


0.44 


0.10 


40( 8) 


3.6 


6.9 


5.3(1.0) 


Chal-C27 


42 


61.1 X 36.7 


8.6x4.5 


1.9 


7.0 


8 


0.042 


0.21 


0.13 


32( 5) 


3.8 


9.1 


4.1(0.7) 


Chal-C28 


43 


59.6 X 32.1 


8.4 X 3.6 


2.3 


6.6 


8 


0.040 


0.17 


0.087 


46( 9) 


3.5 


5.9 


6.0(1.2) 


Chal-C29 


44 


49.3 X 36.1 


6.7 X 4.4 


1.5 


7.1 


9 


0.043 


0.17 


0.13 


33( 6) 


3.8 


8.8 


4.3(0.7) 


Chal-C30 


45 


105.2 X 51.2 


15.5 X 7.0 


2.2 


6.3 


7 


0.038 


0.46 


0.099 


39( 8) 


3.4 


6.7 


5.1(1.0) 


Chal-C31 


46 


71.6 X 44.6 


10.3 X 5.9 


1.7 


6.5 


5.7 


0.039 


0.28 


0.10 


39( 7) 


3.5 


6.9 


5.1(1.0) 


Chal-C32 


47 


91.5 X 49.5 


13.3 X 6.7 


2.0 


6.5 


11 


0.039 


0.40 


0.10 


38( 7) 


3.5 


7.1 


4.9(0.9) 


Chal-C33 


49 


58.2 X 36.0 


8.1 X 4.4 


1.9 


6.5 


8 


0.039 


0.18 


0.10 


38( 7) 


3.5 


7.0 


5.0(1.0) 


Chal-C34 


50 


139.3 X 42.7 


20.7 X 5.6 


3.7 


6.3 


5.5 


0.038 


0.50 


0.10 


37( 7) 


3.4 


6.9 


4.9(1.0) 


Chal-C35 


52 


55.8 X 26.6 


7.7 X 2.4 


3.2 


6.2 


13 


0.038 


0.12 


0.076 


49(11) 


3.4 


5.2 


6.5(1.4) 



rate Ssf in M© yr"' kpc"^. Assuming a mean lifetime in the de- 
tectable prestellar (+ Class 0) stage of 0.5 Myr (see Sect. 17.4b 
and a mean final stellar mass of 0.5 Mq would give a conversion 
factor from 2sf in Mq yr"' kpc"^ to Scores in pc"^ equal to 1. 

To quantify the relationship between star formation and gas 
column density in this local cloud we fit a power law of the form 
Scores = ^2^aas- One has to be very careful with the uncertainty 
budget when dividing two quantities each of which may lie close 
to zero, particularly if this applies for the denominator, and here 
the area of the cloud becomes small at high extinction. We as- 
sume the uncertainties are dominated by Poisson ( V^) counting 
statistics i n each extinction bin and use the Bayesian method as 
laid out in lHatchelL et al.l(l2005l Appendix A) to calculate the re- 
sulting uncertainties on the ratio. The greyscale in Fig.|7j) gives 
the probability distribution for the true value of the mean surface 
density of cores for each Ay value, given the measured counts 
and assuming Poisson statistics. At high extinction values, where 
the number of pixels is close to zero, the mean surface density 
of cores is very uncertain and this shows as a broad spread in the 
greyscale. Using these uncertainties, we perform a minimisation 
to calculate the power-law index a and its associated uncertain- 
ties. Fitting over the full range (A,, = 0-20) gives a steep power- 
law index a - 4.34^Q jg (68.3% confidence limits). This is even 
steeper than the power law found in Perseus (a = 3.0 ± 0.2, 



iHatchell et ai]|2005h . The power-law fit suggests an alternative 
interpretation to that of an Ay threshold, which is that star for- 
mation is simply unusual at low column densities, and a wide 
area search is required to find a dense core wh ich is not em- 
bedded in an extended area of raised extinction (IHatchell et al.l 
2005). 

However, we note that the previous fit tends to underesti- 
mate the mean surface density of cores in the range 8-10 mag 
(see insert of Fig.|7J)), where the peak of the distribution in num- 
ber is located (Fig. |6^). This suggests that the peak at Ay ~8- 
9 mag may have a sp ecific physical signification, as proposed by 
IAndr6etal.l(l2010al) (see Sect.[TJ. 

6.1.2. Sizes 

The source sizes along the major and minor axes before and af- 
ter deconvolution are listed in Cols. 3 and 4 of Table|6] and their 
distribution is shown in Fig. [6^ and b, respectively, along with 
the distribution of mean size (geometrical mean of major and 
minor sizes, i.e. -y/FW//Mmaj x FWHM„^i„). The average ma- 
jor, minor, and mean sizes are 65 '^^^ " , 38 ", and 49 ^]^^ ", re- 
spectively. Only three sources have a major FWHM size larger 
than 110", and no source has a minor or mean FWHM size 
larger than 80". The results of the Monte-Carlo simulations in 
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Table 6. continued. 



Name 




FWHM'' 


FWHM'' 


n h 


^peak 


Ay" 




Mot' 


^50"' 


Cm' 




"50 








(" X ") 


(1000 AU)^ 




(10^' cm" 


'-) (mag) 


(Mo) 


(Mo) 


(Mo) 


(%) 


(10' cm- 


(10-* cm" 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Clial-C36 


53 


49.1 X 33.7 


6.6 X 3.9 


1.7 


6.5 


9 


0.039 


0.14 


0.092 


42( 8) 


3.5 


6.3 


5.5(1.1) 


Clial-C37 


54 


77.6 X 30.1 


11.2 X 3.2 


3.5 


6.1 


7 


0.037 


0.19 


0.091 


41( 8) 


3.3 


6.2 


5.4(1.1) 


Clial-C38 


55 


83.4 X 44.6 


12.1 X 5.9 


2.1 


6.2 


11 


0.038 


0.31 


0.097 


39( 8) 


3.3 


6.6 


5.1(1.0) 


Clial-C39 


56 


82.3 X 35.8 


11.9 X 4.3 


2.8 


6.0 


6 


0.036 


0.24 


0.098 


37( 8) 


3.2 


6.7 


4.8(1.0) 


Clial-C40 


57 


58.7 X 46.6 


8.2 X 6.2 


1.3 


6.1 


12 


0.037 


0.23 


0.095 


39( 8) 


3.3 


6.4 


5.1(1.0) 


Clial-C41 


58 


48.3 X 31.9 


6.5 X 3.6 


1.8 


6.1 


12 


0.037 


0.13 


0.083 


45(10) 


3.3 


5.6 


5.9(1.3) 


Clial-C42 


59 


48.7 X 26.0 


6.6 X 2.3 


2.9 


6.0 


8 


0.036 


0.10 


0.064 


57(14) 


3.2 


4.3 


7.4(1.8) 


Clial-C43 


60 


90.0 X 45.6 


13.1 X 6.0 


2.2 


5.7 


9 


0.035 


0.32 


0.084 


41(9) 


3.1 


5.7 


5.4(1.2) 


Glial -C44 


61 


46.5 X 42.4 


6.2 X 5.5 


1.1 


5.8 


5.9 


0.035 


0.15 


0.074 


48(11) 


3.1 


5.0 


6.3(1.4) 


Chal-C45 


64 


73.2 X 21.4 


10.5 X 2.1 


5.0 


5.7 


5.1 


0.035 


0.12 


0.074 


47(11) 


3.1 


5.0 


6.1(1.4) 


Chal-C46 


65 


36.6 X 21.2 


4.5 X 2.1 


2.1 


5.7 


7 


0.034 


0.059 


0.051 


67(19) 


3.0 


3.5 


8.8(2.4) 


Chal-C47 


67 


87.6 X 36.2 


12.8 X 4.4 


2.9 


5.7 


9 


0.034 


0.24 


0.090 


38( 8) 


3.1 


6.1 


5.0(1.1) 


Chal-C48 


68 


50.7 X 24.0 


6.9 X 2.1 


3.3 


5.5 


7 


0.033 


0.090 


0.059 


56(15) 


2.9 


4.0 


7.4(1.9) 


Chal-C49 


70 


40.2 X 26.9 


5.1 X 2.5 


2.1 


5.6 


9 


0.034 


0.082 


0.072 


47(11) 


3.0 


4.9 


6.2(1.5) 


Chal-C50 


71 


60.6 X 54.3 


8.5 X 7.5 


1.1 


5.5 


12 


0.033 


0.24 


0.086 


39( 9) 


2.9 


5.8 


5.1(1.2) 


Chal-C51 


72 


34.9 X 31.9 


4.2 X 3.6 


1.2 


5.6 


14 


0.034 


0.084 


0.065 


52(13) 


3.0 


4.4 


6.8(1.7) 


Chal-C52 


73 


50.5 X 27.8 


6.9 X 2.7 


2.6 


5.3 


17 


0.032 


0.100 


0.065 


49(13) 


2.9 


4.4 


6.5(1.7) 


Ctiai-Cj J 


/4 


J/.O X 4U. i 


/.I X J. i 




J.J 


t: 
O 


0.033 


0.16 


n no A 
0.0o4 


4U( 9) 




J. / 


D.Z(i.z) 


Clial-C54 


78 


36.6 X 29.9 


4.5 X 3.2 


1.4 


5.2 


17 


0.032 


0.077 


0.060 


53(14) 


2.8 


4.1 


6.9(1.8) 


Clial-C55 


79 


40.7 X 28.8 


5.2 X 2.9 


1.8 


5.3 


10 


0.032 


0.084 


0.077 


42(10) 


2.9 


5.2 


5.5(1.3) 


Clial-C56 


80 


46.3 X 22.4 


6.2x2.1 


2.9 


5.2 


9 


0.031 


0.073 


0.056 


57(16) 


2.8 


3.8 


7.4(2.1) 


Clial-C57 


81 


51.1 X 40.1 


7.0 X 5.1 


1.4 


5.4 


7 


0.032 


0.15 


0.12 


27( 6) 


2.9 


8.1 


3.6(0.8) 


Clial-C58 


82 


48.5x31.6 


6.5 X 3.5 


1.9 


5.5 


11 


0.033 


0.11 


0.088 


38( 9) 


3.0 


6.0 


5.0(1.1) 


Chal-C59 


83 


41.0 X 21.4 


5.3 X 2.1 


2.5 


5.3 


14 


0.032 


0.063 


0.052 


61(17) 


2.9 


3.5 


8.0(2.3) 


Chal-C60 


84 


61.6 X 43.7 


8.7 X 5.7 


1.5 


5.1 


9 


0.031 


0.18 


0.071 


43(11) 


2.8 


4.8 


5.7(1.4) 



Notes. Numbering of Gaussclumps sources like in Table |2] Size of the fitted Gaussian, deconvolved physical source size, and aspect ratio. 
Sizes smaller than 25.4" were set to 25.4" to compute the deconvolved sizes, in order to account for a fit inaccuracy corresponding to a 5cr detection 
in peak flux density. As a result, the minimum size that can be measured is about 2100 AU. The aspect ratio is the ratio of the deconvolved sizes 
along the major and minor axes. '^^ Peak H2 column density computed assuming a dust opacity of 0.01 cm^ g"'. The statistical rms uncertainty is 
1.1 X 10^' cm"-. Visual extinction derived from 2MASS. Mass in the central beam (HPBW = 21.2"), total mass derived from the Gaussian 
fit, and mass computed from the flux measured in an aperture of 50" in diameter. The statistical rms uncertainties of Mpeak and M^^" are 0.006 and 
0.011 Mq, respectively. '•■^^ Mass concentration mpeak/^ijo"- The statistical rms uncertainty is given in parentheses. Beam-averaged free-particle 
density within the central beam and mean free-particle density computed for the mass Mj^// in the aperture of diameter 50". The statistical rms 
uncertainties are 5.6 x 10'' and 7.2 x 10^ cm"^, respectively. Density contrast npcak/^so"- The statistical rms uncertainty is given in parentheses. 



the elliptical case ( Appendix l A. 1 l and Table lA. ft imply that these 
sources, although most of them are faint with a peak flux density 
lower than 150 mJy/beam, are not significantly affected by the 
spatial filtering due to the sky noise removal, with less than 15% 
loss of peak flux density and size. 

The fitting done on the control sample of artificial, circular, 
Gaussian sources shows that the accuracy of the sizes derived 
for sources with a peak flux equal to ~ 5cr is roughly 20% (see 
Figs. lA. It and d). With HPBW - 21.2", the size accuracy of a 
weak ~ 5a- unresolved source is 4.2". Therefore, faint sources 
with a size smaller than ~ 25.4" cannot be reliably deconvolved 
and we artificially set their size to 25.4" to perform the deconvo- 
lution. As a result the minimum deconvolved FWHM size that 
we can measure is ~ 14" (2100 AU). The average deconvolved 
mean FWHM size is 43 ", i.e. 6500 AU (see Fig. 

lEnoch et"an (|2008) measure average deconvolved mean 
sizes of 16500 + 6000, 10700 ±3100, and 7500 ± 3500 AU 
(FWHM) for the starless cores extracted with Clumpfind from 
their Bolocam maps of Perseus, Serpens, and Ophiuchus, re- 
spectively (see Col. 9 of Table Q. With the MPIfR 19-channel 
bolometer array of the IRAM 30m telescope {HPBW - 15" af- 
ter smoothing), 'Motte et aP (1199 8) derive significantly smaller 
sizes for the starless cores in Ophiuchus based on a method sim- 
ilar to the one used here for Cha I (multiresolution decompo- 
sition and Gaussian fitting): they find an average deconvolved 



mean size of 2100+}^gg AU when rescaled to the same distance 
as assumed bv lEnoch et alj (12008). Both the differences in angu- 
lar resolution (31" versus 15") and in the extraction method used 
(no background subtraction versus multiresolution analysis) cer- 
tainly explain the discrepancy of a factor of 3 between both 
studies. Using Gaussclumps without background subtraction, 
Curtis & Richer (2010) derive an average deconvolved mean 
size of 3800 + 200 AU for the starless cores observed in Perseus 
by Hatchell et al.. (.2005.) with SCUBA, i.e. a factor of 4 times 
smaller than the Bolocam cores. The angular resolution may 
therefore be the dominant factor influencing the derived FWHM 
sizes. The angular resolution of LABOCA (21" after smooth- 
ing) being a bit larger than the SCUBA and 30 m resolutions, we 
can conclude that the Cha I starless cores have similar physical 
sizes than the Perseus cores, are probably larger than the Serpens 
cores (maybe by a factor of 1 .5-2), and are certainly larger than 
the Ophiuchus cores (by a factor of 2-3). 

No complete survey of Tau rus in dust continuum emission 
exists. iKauffmann et al.l (l2008h mapped 10 dense cores with 
MAMBO at the IRAM 30 m telescope. After smoothing to 20", 
they identified 23 starless subcores. The average deconvolved 
mean FWHM size of this sample is 20600+;^^°° AU (see TableHl. 
If this small sample is representative for the entire population of 
starless cores in Taurus, then the Cha I starless cores are signifi- 
cantly smaller (by a factor of 3) than the Taurus ones. 
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Fig. 6. Distribution of physical properties obtained for the 60 starless (or Class 0) sources found with Gaussclumps in the sum map 
of Cha 1 at scale 5. The mean, standard deviation, and median of the distribution are given in each panel. The asymmetric standard 
deviation defines the range containing 68% of the sample, a FWHM sizes along the major (solid line) and minor (dashed line) axes. 
The filled histogram shows the distribution of geometrical mean of major and minor sizes. The mean and median values refer to the 
filled histogram. The dotted line indicates the angular resolution (21.2"). b Same as a but for the deconvolved sizes, c Aspect ratios 
computed with the deconvolved sizes. The dotted line at 1.4 shows the threshold above which the deviation from 1 (elongation) can 
be considered as significant, d Peak H2 column density. The dotted line at 5.25 x 10^' cm"^ is the 5cr sensitivity limit, e Visual 
extinction derived from 2MASS. 



6.1.3. Aspect ratios 

The distribution of aspect ratios computed with the deconvolved 
FWHM sizes is shown in Fig. |6};. Based on the Monte Carlo 
simulations of Appendix lA.il we estimate that a faint source can 
reliably be considered as intrinsically elongated when its aspect 
ratio is higher than 1 .4. 17% of the sources are below this thresh- 
old while 83% can be considered as elongated. The average as- 
pect ratio is 2.1 ^g g. It is similar to the ones measured in Serpens 
and Taurus, somewhat larger than in Perseus, and significantly 
larger than in Ophiuchus (see Col. 10 of Table |7]l. 



6.1.4. Column densities 



The average peak H2 column density of the starless sources in 
Cha I is 9 +4 X 10^' cm"^ (Fig.|6}l). This is 4 times lower than the 
average peak column density of the starless cores in Perseus and 
Serpens, and 3 to 7 times lower than in Ophiuchus (see Col. 1 1 
of Table |7] and Appendix IB. 61 for details). It appears to be sig- 
nificantly smaller than in Taurus too (by a factor of 2.6), but 
since the Taurus sample is not complete and the source extrac- 
tion methods differ, this may not be significant. 



6.1 .5. Masses and densities 

The distribution of masses and free-particle densities listed in 
Table |6] are displayed in Fig.|8] The 5cr sensitivity limit used to 
extract sources with Gaussclumps corresponds to a peak mass of 
0.032 M0 and a peak density of 2.8 x 10^ cm""*, computed for 
a diameter of 21.2" (3200 AU). The median of the peak mass 
distribution is 0.039 M©, implying a median peak density of 
3.5 X 10'' cm"^ (Figs. [8^ and e). We give the mass integrated 
within an aperture of diameter 50" (7500 AU) in Col. 10 of 
Table |6] This aperture is well adapted to our sample for three 
reasons: it corresponds to the average mean, undeconvolved 
FWHM size (see Sect. l6.1.2l and Fig.|3i), it is not affected by the 
spatial filtering due to the sky noise removal (see Appendix lA.il 
and Table lA.U . and it is still preserved in the sum map at scale 
5 (see Appendix ICl and Table ICTT i. The median of the mass in- 
tegrated within this aperture is 0.10 M©, corresponding to a me- 
dian mean density of 7 x 10"^ cm"^ (Figs. [SJ) and f). 

The cores in Cha 1 are about 3^ times less dense within 
an aperture of diameter 7500 AU than those in Perseus and 
Serpens, 6 times less dense than those in Ophiuchus, and 3 
times less dense than those in Taurus (see Col. 12 of Table [T] 
and Append ix IB. 6. 21 f or details). H owever, both studies of 
Eno ch et all 12008) and lKauffmann et al. (2008) do not filter the 
extended emission as we do with our multiresolution decompo- 
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Fig. 7. a Distribution of extinction across the map. The area for Ay < 3 regions is Hmited by the map boundary, b Mean surface 
density of starless (+ Class 0) cores as a function of gas column density as measured by Ay- For each Ay value, the probability 
distribution for the true value of the mean surface density is displayed in greyscale. The corresponding colorbar on the right is 
labeled in natural logarithm. The insert is a zoom on the low-extinction range. The best-fit power-law curve to the mean surface 
density is shown in green/black (see Sect. l6.1.Tl for details). 



Table 7. Sensitivities of (sub)mm surveys of nearby molecular clouds and average properties of the detected starless cores. 



Cloud 

(1) 


Dist. 

(pc) 

(2) 


Freq. 

(GHz) 

(3) 


HPBW 

(") 
(4) 


{mJy/beam) 

(5) 


o-n' 

(10-' cm"- 

(6) 


J500 c 
" M 
) (Mo) 

(7) 


Ref." 

(8) 


< FWHM >' 
(1000 AU) 

(9) 


< Ra >■' 
(10) 


< A'peak 
(10-' cm"-) 

(11) 


< "7500 AU >'' 
(10' cm"^) 

(12) 


Ref." 
(13) 


Chamaeleon I 
Ophiuchus 


150 
125 


345 
268 


21 
31 


12 
27 


1.1 
2.4 


0.012 
0.035 


1 

2 


6 5+3° 
7.5 ±3.5 


9 1+0.8 
^•^-0.6 

1.3 ±0.2 


9+5 

60 


1 0+'-' 
5.7-6.5 


1 

8 






240 


13 


8 


5.9 


0.035 


3 




1 3+°-^ 

'^■-'-0.3 


26+23 

^"-13 




3 


Taurus 


140 


250 


11 


1.4 


1.3 


0.007 


4 


20.6!« 


1 9+0.5 


23+13 


3.0 


4 


Perseus 


250 


268 


31 


15 


1.4 


0.039 


5 


16.5 ± 6.0 


1.7 ±0.9 


38 


3.5-3.9 


8 






353 


14 


35 


6.6 


0.086 


6 


3.8 ±0.2 


1.4 ±0.1 


35 ±3 




9 


Serpens 


260 


268 


31 


9.5 


0.9 


0.026 


7 


10.7 ±3.1 


1.7 ±0.5 


33 


2.8-3.2 


8 



Notes. The values in Cols. 6 and 7 were computed with the same assumptions for all clouds (see Appendix iBt. The values in Cols. 9 to 12 
obtained from the litterature were rescaled to match the physical assumptions and size criteria adopted for our analysis of Cha 1 (see Appendix lB.6l 
for details). Peak flux density rms sensitivity in the original beam. Peak column density rms sensitivity in the original beam. Mass rms 
sensitivity per aperture of diameter 7500 AU. References: (1) this work: (2)[Youn g et al. (2006); (3) Motte et al. ( 1998); (4) Kau ffmann et al.l 
(l2008h : (5) lEnochet an (l2006): f6 VHatchell et al.l ( l2005l) : (7) lEnoch et aP ( 120071) : (8) Enoch et al. (2008); (9) Curtis & Richer (2010). Average 
deconvolved mean size of the detected starless cores. Average aspect ratio. Average peak column density in a FWHM size of 3180 AU, 
corresponding to our angular resolution of 21.2" at the distance of Cha I. Average mean density in an aperture of diameter 7500 AU. 



sition to isolate the starless cores from the environment in which 
they are embedded. This may bias the densities in both stud- 
ies to somewhat higher values compared to Cha I. In addition, 
the sensitivities of the different surveys in terms of mass (and 
mean density) within the aperture of 7500 AU are not the same 
(see Col. 7 of Table |7]i. Our survey is more sensitive than those 
reported by Enoch et a l. (2008) by a factor of 2 to 3, but less 
sensitive than the one of Kau ffman n et al.l (120081) bv a factor of 
nearly 2. This means that the surveys of Perseus, Serpens, and 
Ophiuchus may be biased towards higher masses and mean den- 
sities and that they may have missed a population of less dense 
starless cores similar to the one found in Cha I. In the case of 
Taurus, the bias of the source sample of iKauffmann et all (l2008l) 
prevents any conclusion concerning the possible existence of a 
population of less dense starless cores. 



6.3cr for the average size of the source sample (FWHM = 49"fl 
It coiTesponds to a total mass of 0.22 Mq. The median total mass 
is very close (0.21 Mq), which implies that only 50% of the de- 
tected sources are above the estimated 90% completeness limit. 
For comparison, the 90% completeness limit of the starless core 
sample of Enoch et al. (2008) in Perseus is ~ 0.9 Mq, i.e. 1.2 Mq 
when rescaled to the same temperature and dust opacity as we 
use here (12 K and /3 = 1.85). Our complete ness limit is sim- 
ilar to that obtained by iKonvves et alj (1201 Oh for their 1 1 deg^ 
sensitive continuum survey of the Aquila Rift cloud complex 
(distance 260 pc) with Herschel: assuming a mass distribution 
proportional to the radius, they estimated their sample of 541 
starless cores to be 75 and 85% complete above masses of 0.2 
and 0.3 M©, respectively. Their assumed dust opacity law yields 
""870 - 0.012 cm^ g"', i.e. only 20% different from our assump- 
tion so the completeness limits can be directly compared. 



Figure |8}; shows the distribution of total masses computed 
from the Gaussian fits (Col. 9 of Table |6]l. The completeness 
limit at 90% is estimated from a peak flux detection threshold at 



* For a Gaussian distribution of mean value m and standard deviation 
cr, the relative population below m - 1.28cr represents 10%. Therefore 
our peak flux detection threshold at 5cr implies a 90% completeness 
limit at 6.3cr, with cr the rms noise level in the sum map at scale 5. 
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P 

1 1.5 2 2.2 2.4 2.6 




5.5 6 6.5 4.5 5 5.5 6 4 6 8 10 

logio(npeak (cm ^)) logio(n5o" (cm ^)) n^^^y-/n^g.. 

Fig. 8. Distribution of masses and densities obtained for the 60 starless (or Class 0) sources found with Gaussclumps in the sum map 
of Cha I at scale 5. The mean, standard deviation, and median of the distribution are given in each panel. The asymmetric standard 
deviation defines the range containing 68% of the sample, a Peak mass of the fitted Gaussian, b Mass within an aperture of diameter 
50". c Total mass of the fitted Gaussian. The dotted line indicates the estimated completeness limit at 90% for Gaussian sources 
corresponding to a 6.3cr peak detection limit for the average source size, d Mass concentration, ratio of the peak mass to the mass 
within an aperture of diameter 50". e Peak density, f Mean density within an aperture of diameter 50". g Density contrast, ratio of 
peak density to mean density. In panels a and e, the dotted line indicates the 5cr sensitivity limit. The upper axis of panel d, which 
can also be used for panel g, shows the power-law exponent p derived assuming that the sources have a density profile proportional 
to r and a uniform dust temperature. Alternately, the upper axis of panel g, which can also be used for panel d, deals with the 
case where the density is uniform within a diameter Dflat and proportional to r"^ outside, still with the assumption of a uniform 
temperature. See Sect. l6.1.5] for the Umitations of these upper axes. 



We estimate the mass concentration of the sources from the 
ratio of the peak mass to the mass within an aperture of 50" 
(Col. 11 of Table |6]l which is insensitive to the spatial filtering 
due to the data reduction (see Table IaTTT i. A similar property is 
the density contrast measured as the ratio of the peak density 
to the mean density within this aperture (Col. 14 of Table |6]l. 
The statistical rms uncertainties on the peak mass and the mass 
within 50" are 0.006 and 0.011 M©, respectively, which means 
a relative uncertainty of up to 20% for the weakest sources. The 
distributions of both ratios are shown in Figs.[8}l and g and their 
rms uncertainties are given in parentheses in Cols. 1 1 and 14 
of Table |6] The two outliers with the largest ratios are also those 
with the highest relative uncertainty (about 30%). The upper axis 
of Fig.[8}l, which can also be used for Fig. [8^, displays the expo- 
nent of the density profile under the assumptions that the sources 
are spherically symmetric with a power-law density profile, i.e. 
p cx r^P, and that the dust temperature is uniform. The median 
mass concentration and density contrast are 0.39 and 5.1, respec- 

' The relative uncertainty of the ratio is equal to the square root of 
the quadratic sum of the relative uncertainties of its two terms, i.e. we 
assume both terms are uncorrelated. 



tively. This corresponds to p ~ 1 .9, suggesting that most sources 
are significantly centrally-peaked. It is very close to the expo- 
nent of the singular isothermal sphere {p - 2). As a caveat, we 
should mention that the relation used to derive p is valid only 
for sources with a mass between 21.2" and 50" that is mea- 
surable and extends up to 50". In practice, most sources have 
Cm < 50% with a relative uncertainty less than 20% (see Col. 1 1 
of Table |6]l, meaning that for those sources, the mass between 
both diameters is measured with a reasonable accuracy. The few 
sources with Cm > 50% (13% of the sample), especially those 
with Cm > 60%, are also those with the highest relative uncer- 
tainties (~ 30% for the latter). The relation to derive p may not 
be valid for these few sources, or they are simply too weak to 
give a reliable estimate. As a result, these sources may bias the 
mean and median values of p to slightly higher values, but since 
they are not numerous, the bias is certainly smal0. 

The upper axis of Fig. [8^, which can also be used for Fig.[8}l, 
deals with an alternate case where the sources have a constant 



If we remove the sources with Cm > 60% and Cm > 50%, the mean 
of Cm is only reduced by 2% and 6%, respectively, and the median 
decreases by less than 1%. 
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density within a diameter Dgat and a density proportional to 
r"^ outside, still with the assumption of a uniform temperature. 
Under these assumptions, the measurements are consistent with 
a flat inner region of diameter 16" at most (2400 AU) for a few 
sources, but most sources have Dgat < 10" (1500 AU), or cannot 
be described with such a density profile. 



6.1 .6. Mass versus size 

The distribution of total masses versus source sizes derived 
from the Gaussian fits is shown in Fig. |9^. About 50% of the 
sources are located between the Scr detection limit (full line) 
and the estimated 90% completeness limit (dashed line), sug- 
gesting that we most likely miss a significant number of sources 
with a low peak column density. Figure |9j5 shows a similar di- 
agram for the deconvolved source size. If we assume that the 
deconvolved FWHM size is a good estimate of the external 
radius of each source, then we can compare this distribution 
to the critical Bonnor-Ebert mass that characterizes the limit 
above which the hydrostatic equilibrium of an isothermal sphere 
with thermal support only is gravitationally unstable. This re- 
lation Mbe(/?) = 2ARallG (lBonnodll95"6i) . with Mbe(^) the 
total mass, R the external radius, as the sound speed, and G the 
gravitational constant, is drawn for a temperature of 12 K as a 
solid line in Fig. Only three sources (Chal-Cl to C3) are 
located above this critical mass limit. If we account for a fac- 
tor of 2 uncertainty on the mass, then seven additional sources 
could also fall above this limit provided the mass is underesti- 
mated (Chal-C4 to C7, C9, Cll, and C13, see dash-dotted line 
in Fig. Ig]?). Most sources, however, have a mass lower than the 
critical Bonnor-Ebert mass by a factor of 2 to 5. The uncertainty 
on the temperature (see Appendix IB. 2l i does not influence these 
results much since, even in the unlikely case of the bulk of the 
mass being at a temperature of 7 K, the measured masses would 
move upwards relative to the critical Bonnor-Ebert mass limit by 
a factor of 1 .9 only, because the latter is also temperature depen- 
dent. 

Following the analysis of iKonvves et al.l (1201 Ol) for the 
Aquila starless cores, we can also estimate the critical Bonnor- 

Ebert mass w ith the equation MBE(^'ext) = l-l^ '^s ^ext 
(lBonnoill956l) . The external pressure is estimated with the equa- 
tion Pext - 0.88 G (fimiiNci)^, with /i the mean molecular weight 
per free particle, mu the mass of hydrogen, and A^ci the column 
density of the local ambient cloud in which the sources are em- 
bedded (.McKee & Tan ,2003). We use the extinction listed in 
Table|6]to estimate this background column density with the con- 
version factors given in Appendix IB. 3 1 The sources with a mass 
larger than MBE(^'ext) are the same as for Mbe(^), plus Chal-C4. 
The agreement between both estimates of Mbe suggests that our 
estimates of the external radius and external pressure are consis- 
tent. If we account again for a possible factor of 2 uncertainty on 
the mass measurement, only three additional sources would fall 
above the critical limit based on the pressure (Chal-C7, Cll, 
and CI 8). In summary, four sources are likely above the criti- 
cal Bonnor-Ebert limit (Chal-Cl to C4), and seven additional 
ones may also be if their mass is underestimated by a factor of 
2 (Chal-C5 to C7, C9, CI 1, C13, and C18). The implications of 
this analysis will be discussed in Sect.Q 

The mass concentration Cm is plotted versus source size in 
Fig Et. Cm is actually equal to the ratio of the peak flux to the 
flux integrated within the aperture of diameter 50". When the 
sources do not overlap, this ratio is nearly independent of the 
Gaussian fitting since the second and third stiffness parameters 
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Fig. 10. a Total mass versus visual extinction Ay for the 60 star- 
less (or Class 0) sources found with Gaussclumps in the sum 
map of Cha I at scale 5. The dashed line shows the estimated 
90% completeness limit (0.22 Mq). b FWHM size versus visual 
extinction. The angular resolution (21.2") is marked by the dot- 
ted line. 



of Gaussclumps were set to 1, i.e. biasing Gaussclumps to keep 
the fitted peak amplitude close to the observed one and the fitted 
center position close to the position of the observed peak. The 
dashed line shows the expected ratio if the (not deconvolved) 
sources were exactly Gaussian and circular and allows to esti- 
mate the departure of the sources from being Gaussian within 
50". Most sources have a mass concentration consistent with 
the Gaussian expectation, but many of them have a significant 
uncertainty on Cm that prevents a more accurate analysis. The 
two obvious outliers toward the lower left are sources Chal-C6 
and C8, which have strong neighbors significantly contaminat- 
ing their flux within 50" (sources Chal-Sl and S4, respectively). 

There is no obvious correlation between the total mass or 
FWHM size of the sources and the visual extinction of the en- 
vironment in which they are embedded (see Fig. [TOt . A simi- 
lar conclusion was drawn by ,Sadavov et al. (2010) for the five 
nearby molecular clouds Ophiuchus, Taurus, Perseus, Serpens, 
and Orion based on SCUBA data. 



6.1.7. Core mass distribution (CIVID) 

The mass distribution of the 60 starless (or Class 0) sources is 
shown in Fig.[TT] in both forms AN I AM (a) and AN I A log(M) (b). 
Its shape in Fig. [TTk looks very similar to the shape of the mass 
distribution found in other star forming regions with a power- 
law-like behavior at the high-mass end and a flattening toward 
the low-mass end. In our case, the flattening occurs below the 
estimated 90% completeness limit (0.22 Mq) and may not be 
significant. Above this limit, the distribution is consistent with a 
power-law. The exponent of the best power-law fit {a - -2.29 ± 
0.30 for AN I AM, a\o 



-1.2 9 + 0. 30 for dA^/dlog(M)) is very 



close to the value o f lSalDeted(ll955h that characterizes the high- 
mass end of the stellar initial mass function {a - -2.35) and 
steeper than the expon ent of the typical mass spectrum of CO 
clumps {a = -1.6, see lBlitzl[T993[ iKramer et alJll998h . Such a 
result was also obtained fo r the population o f starless cores in, 

e.g., Ophiuchus (a 2.5, Motte et al.l ll998l) . Serp ens, Perseus 

(a = -2.3 + 0.4. lEnoch et al.' 2008). the Pipe nebula (lAlves et al.1 
'20071), Taurus (ciog = -1.2 + O.fl Sadavov et al."2010), and flie 
Aquila Rift (ffiog = -1.5±0.2. lAndre et al. 2010a: Konvves et al.l 
120101) . 



' Given the small size of the Taurus sample (69 starless sources), the 
uncertainty seems underestimated compared to the other clouds. 
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Fig. 9. a Total mass versus mean FWHM size for the 60 starless (or Class 0) sources found with Gaussclumps in the sum map 
of Cha I at scale 5. The angular resolution (21.2") is marked by the dotted line. The solid line (M oc FWHM^) is the 5cr peak 
sensitivity limit for Gaussian sources. The dashed line shows the 6.3cr peak sensitivity limit which corresponds to a completeness 
limit of 90% for Gaussian sources, b Total mass versus mean deconvolved FWHM size. Sizes smaller than 25.4" were set to 
25.4" before deconvolution (see note b of Table|6l). The solid line shows the relation M - lARa^/G that characterizes critical 
Bonnor-Ebert spheres (see Sect. 16. 1.6b . The dash-dotted and dashed lines show the location of this relation when divided by 2 and 5, 
respectively. The sources with a mass larger than the critical Bonnor-Ebert mass Mbe.p^i estimated from the ambient cloud pressure 
are shown with a filled circle. Those with 0.5 Mbea^ < M < Mbe,p„, are shown with a circle, c Mass concentration versus mean 
FWHM size. The dashed line is the expectation for a circular Gaussian flux density distribution. 



The mass distribution of the Cha I starless sources does 
not show any significant flattening down to the estimated 90% 
completeness limit (0.22 Mq), although the signal-to-noise ra- 
tio may not be sufficient to draw a firm conclusion. This 
seems to contrast with the mass distributions of Ophiuchus and 
Aquila which flatten below 0.4 and 1.0 Mq, a factor ~ 4 and 
~ 3 above their completene ss limits, respectively (Mot te et al.l 
119981: iKonvveset all 1201 Ol) . This flattening translates into a 
dA^/dlog(M) curve with a significant maximum around 0.5- 
0.6 Mq for Aquila while, in Ophiuchus, it keeps increasing with 
a reduced expone nt (a\og ~ -0.5) down to the completeness 
Hmit (~ 0.1 Mo, lMotte etal.|[l998l) . The dA^/dlog(M) curve 
in Cha I reaches a maximum at a mass of 0.2 M p (Fig. fTTb). 
significantly lower than in Aquila. ISadavov et al.l (POlO) found 
a maximum around 0.16 Mq for Taurus, which a priori looks 
similar to our result for Cha I. However, their analysis is based 
on masses computed within a 850 fim contour of 90 mJy/14"- 
beam, meaning that they miss a large fraction of the total mass 
of each core. It is very likely that the mass corresponding to 
the true maximum is higher than 0.16 Mq by a factor of a few, 
even maybe an or der of magnitude. For comparison, the average 
mass measured bv lKauffmann et all (|2008) within an aperture of 
diameter 8400 AU, which represents only a fraction of the to- 
tal mass (roughly 50%), is already a factor of 3 higher for their 
sample of 28 starless "peaks" (see Sect. 16.1. St . The discre pancy 
with the mass distribution obtained bv .Onishi et'aTI (|2002|) based 
on H'^CO^ observations is even larger, since their distribution 
dN/dM flattens around 2 Mq. Therefore, the true maximum of 
the dA^/dlog(M) curve in Taurus is certainly at a mass signifi- 
cantly higher than in Cha I. 

6.2. YSO candidates 

The sources extracted with Gaussclumps in the sum map at 
scale 5 and associated with known YSOs found in the SIMBAD 
database are listed in Table [8] Their possible SIMBAD asso- 
ciations are listed in Table |2] Most of these 16 sources are 
barely resolved. Their deconvolved sizes are computed from the 
fitted sizes after multiplying the latter by (1 -i- l/SNR), with 



SNR the peak signal-to-noise ratio. In this way, the deconvolved 
sizes are about Icr upper limits on the true physical sizes of 
the barely resolved sources. The peak masses, total masses, and 
masses within an aperture of diameter 50" are computed with 
the same assumptions as for the starless sources, except for the 
dust temperature and mass opacity that we here take as 20 K and 
0.03 cm^ g"', respectively (see AppendixlBl). The mass concen- 
tration Cm is higher than 70% for 1 1 sources (69%), consistent 
with those sources being nearly point-like and dominating their 
environment within a diameter of 50". The remaining 5 sources 
are contaminated by nearby sources (case of Chal-S2 and S4) 
or were qualified as "candidate" (Sc) in Table |2] because they 
are significantly embedded in a larger-scale dense core and are 
difficult to extract from their environment (Chal-S13 to S15). 

The same analysis was performed for the 5 additional com- 
pact sources found in Sect. l5.4l based on a SIMBAD association. 
Their deconvolved sizes, peak column density, and masses are 
given in Table |9] They all have a concentration parameter con- 
sistent with their being nearly point-like, but the uncertainties 
are large because they are faint and close to the detection limit. 

In total, 21 compact sources associated with a YSO in 
the SIMBAD database are detected above the 3.5o- level (42 
mJy/21.2"-beam) in the 870 fim continuum map of Cha I. This 
represents about 10% of t he known YSO members of Cha I (237 
YSOs, see lLuhmanl2008 h. The infrared class of the nearest YSO 
associated to each compact source is given in Tables [8] and |9] 
There are 3 Class I, 2 "Flat", and 16 Class II sources. This repre- 
sents a 870 pm detection rate of 75%, 20%, and 17% for the re- 
spective infrared classes, based on the census o f lLuhmanl(l2008h 
(see Table [TTli. Two of the detected Class II sources are transi- 
tional disks, i.e. objects for which the spectral energy distribu- 
tion indicates centr al clearings or gaps in the dust distribution 
(Cha-S5 and S6, see lKim et al] l2b09). The missing Cla ss I YSO 
not listed in Tables[8]and|9](2M Jl 1092855-7633281 in|LuhmaB 
2008) is located about 8" east of Cha-C12. Given that Cha-C12 
is extended and does not show any peak at the position of the 
YSO, we did not consider the 870 pm emission as being primar- 
ily associated with the YSO itself. We can however not exclude 
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Table 8. Characteristics of YSOs extracted with Gaussclumps in the 870 fim continuum map of Cha I filtered up to scale 5. 



Name 




FWHM" 


FWHM" 


SNR' 


N peak 










Cm' 


Class" 


Ref." 




(" X ") 


(1000 AU)2 




(10^' cm-2) 


(mag) 


(Mo) 


(Mo) 


(Mo) 


(%) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Chal-Sl 


1 


23.5 X 22.9 


1.6 X 1.3 


215 


33 


7 


0.20 


0.24 


0.23 


86( 1) 


II 


1 


Chal-S2 


2 


26.6 X 26.6 


2.5 X 2.5 


101 


15 


9 


0.092 


0.14 


0.13 


69( 1) 


I 


1 


Cha 1 -S3 


3 


25.1 X 23.0 


2.1 X 1.4 


97 


15 


11 


0.089 


0.11 


0.11 


79( 1) 


II 


1 


Chal-S4 


5 


57.8 X 31.0 


8.3 X 3.5 


51 


1.1 


14 


0.046 


0.18 


0.097 


48( 1) 


I 


1 


Chal-S5 


7 


22.7 X 21.3 


1.5 X 1.0 


24 


3.1 


2.7 


0.022 


0.024 


0.024 


95(7) 


II 


2 


Chal-S6 


10 


21.4 X 21.2 


1.2 X 1.1 


16 


2.5 


11 


0.015 


0.015 


0.013 


112(15) 


II 


1 


Chal-S7 


12 


21.9 X 21.2 


1.5 X 1.2 


14 


2.1 


6 


0.013 


0.013 


0.013 


99(14) 


Flat 


1 


Chal-S8 


17 


22.0 X 21.2 


1.7 X 1.4 


11 


1.7 


3.1 


0.010 


0.011 


0.0098 


107(19) 


II 


3 


Chal-S9 


18 


21.8 X 21.2 


1.6 X 1.4 


11 


1.6 


7 


0.0098 


0.010 


0.0095 


103(19) 


II 


1 


Chal-SlO 


20 


21.6 X 21.2 


1.6 X 1.5 


9.9 


1.5 


14 


0.0091 


0.0093 


0.012 


73(12) 


Flat 


1 


Chal-Sl 1 


22 


21.4 X 21.2 


1.6 X 1.5 


9.5 


1.4 


3.2 


0.0087 


0.0088 


0.0078 


112(25) 


I 


1 


Chal-S12 


28 


21.4 X 21.2 


1.8 X 1.7 


1.1 


1.2 


8 


0.0071 


0.0072 


0.0071 


99(25) 


II 


1 


Chal-S13 


31 


35.8 X 28.5 


5.2x3.6 


1.1 


1.2 


7 


0.0071 


0.016 


0.014 


52( 9) 


II 


1 


Chal-S14 


62 


21.5 X 21.2 


2.0 X 1.9 


5.9 


0.89 


11 


0.0054 


0.0054 


0.020 


27( 5) 


II 


1 


Chal-S15 


63 


54.5 X 24.3 


9.1 X 2.9 


5.6 


0.85 


8 


0.0052 


0.015 


0.0097 


53(13) 


II 


4 


Chal-S16 


77 


24.8 X 21.2 


3.1 X 2.1 


5.1 


0.78 


6 


0.0047 


0.0055 


0.0053 


89(31) 


II 


1 



Notes. Numbering of Gawssdumps sources like in Table[2] Size of the fitted Gaussian and deconvolved source size. The deconvolved size 
is computed from the fitted size multiplied by (1 -I- ijSNR), with SNR the peak signal-to-noise ratio, and therefore is an upper limit in most 
cases. Peak signal-to-noise ratio of the fitted Gaussian. Peak H2 column density computed assuming a dust temperature of 20 K and a dust 
opacity of 0.03 cm^ g"'. The statistical rms uncertainty is 1.5 x 10^" cm~^. Visual extinction derived from 2MASS. Mass in the central beam 
(HPBW = 21.2"), total mass derived from the Gaussian fit, and mass computed from the flux measured in an aperture of 50" in diameter. The 
statistical rms uncertainties of Mpeak and M^^i' are 0.0009 and 0.0015 Mo, respectively. Mass concentration mpeak /^i^n". The sta tistical rms 
uncertainty i s given in parenth eses. Infrare d cl ass of the nearest Y SO listed in Col. 11 of Tableland reference: (l) lLuhman et al.l j2008h : (2) 
iFurlan et al.l ( [2009,) : (3) Henning et alj ( 119931) : (4) lLuhman & Muench ( .2008.) . 



Table 9. Characteristics of the additional compact sources listed in Table |5] 



Name 


A'add" 


FWHM" 


FWHM" 


SNR" 


^ peak 


A/ 


Mpeak-' 




M >i > 
'"50 


Cm' 


Class" 


Ref." 






(" X ") 


(1000 AU)- 




(10-' cm-2 


) (mag) 


(Mo) 


(Mo) 


(Mo) 


(%) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Chal-S17 


1 


24.4 X 20.7 


3.1 X 2.1 


4.6 


0.69 


2.7 


0.0042 


0.0047 


0.0048 


87(34) 


II 


2 


Chal-S18 


2 


56.9 X 18.6 


11.3 X 2.1 


2.7 


0.40 


2.3 


0.0024 


0.0057 


0.0039 


63(35) 


II 


2 


Chal-S19 


3 


23.2 X 18.4 


3.0 X 1.4 


3.9 


0.59 


6 


0.0036 


0.0034 


0.0030 


118(66) 


II 


1 


Chal-S20 


4 


20.4 X 16.6 


1.9 X 0.0 


4.9 


0.74 


2.9 


0.0045 


0.0034 


0.0030 


149(82) 


II 


1 


Chal-S21 


5 


18.6 X 16.3 


1.0 X 0.0 


5.1 


0.77 


2.2 


0.0047 


0.0031 


0.0022 


211(152) 


II 


1 



Notes. Numbering of additional sources like in Table[5] Size of the fitted Gaussian and deconvolved source size. The deconvolved size is 
computed from the fitted size multiplied by (1 + I /SNR), with 5A'^ the peak signal-to-noise ratio, and therefore is an upper limit in most cases. 
A zero value means that the uncertainty was still underestimated. Peak signal-to-noise ratio of the fitted Gaussian. '''' Peak H2 column density 
computed assuming a dust temperature of 20 K and a dust opacity of 0.03 cm- g"'. The statistical rms uncertainty is 1.5 x 10-" cm"-. Visual 
extinction derived from 2MASS. Mass in the central beam (HPBW = 21.2"), total mass derived from the Gaussian fit, and mass computed from 
the flux measured in an aperture of 50" in diameter. The statistical rms uncertainties of Mpeak and Mj^/' are 0.0009 and 0.0015 Mo, respectively. 
*'^' Mass concentration m^ cak/tn^n" ■ The statistic al r ms uncertainty is giv en in parentheses. Infrared class of the nearest YSO listed in Col. 10 of 
Table|5]and reference: (D iLuhman eTH] ( 120081) : (2) lFurlan etal] ( 120091) . 



that part of this emission comes from the disk and/or the residual 
envelope of this YSO. 

6.3. Large scale structures 

The smoothed map at scale 5 clearly shows the two main fila- 
mentary structures in fields Cha-Center and Cha-North akeady 
mentioned in Sect. 14. II (see Fig. [Sj)). It also reveals two fainter 
filaments that seem to be connected with a T-shape in the south- 
western part of field Cha-Center. Two other faint, elongated, 
large-scale structures are present in fields Cha-Center and Cha- 
South, respectively. These filamentary structures were akeady 
noticed by Haikala et al. (2005) in their C'^0 1-0 map where 
they appear to be broader, certainly because their map is more 
sensitive to lower densities than our 870 fim dust emission map. 



These authors mentioned that the C'^O filaments are oriented at 
roughly right angles to the magnetic field (see their Fig. 8). 

The properties of the filamentary structures detected with 
LABOCA were obtained with Clumpfind. Since we are here 
interested in the largest structures, only one contour level at 
10 mJy/21"-beam (about 3.3cr in the smoothed map) was de- 
fined to run Clumpfind. In this way, these filaments were not 
splitted into sub-structures. Clumpfind extracted 1 1 structures. 
Only those with a peak flux density higher than 40 mJy/beam 
were selected. Their properties are listed in Table [10] and they 
are labeled in Fig. [S^ accordingly. The major and minor sizes 
given in Table [TO] correspond to the full extent of each filament 
measured by hand using the contour at 10 mJy/beam. 

The length of the filaments varies between 6' and 22' (0.25 
and 1 pc) and their width between 2.3' and 4.9' (0.1 and 0.2 
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Table 10. Large-scale filamentary structures in Cha I probed with LABOCA at 870 /jm (see Fig. |5}'). 
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Notes. Peak position and flux density measured with Clumpfind. Mean and total flux densities inside the contour at 10 mjy/21"-beam (about 
3.3cr) - hereafter called contour C. Angular area of contour C. Major and minor angular axes of contour C measured by hand. Aspect 
ratio. Peak and mean H2 column density. Mass within contour C in smoothed map at scale 5. Physical major and minor axes of contour C. 

Mass within contour C in full LABOCA map. Mass per unit length M,o(/a„,iij. **' This structure has a T shape and consists of two filaments that 
were identified as a single structure by Clumpfind. The major and minor axes were measured for each filament separately, but all other parameters 
correspond to the entire structure. 



pc). Their aspect ratio ranges from 2 to 7. The Monte Carlo sim- 
ulations performed in Appendix I A. II (see Table IATI ) show that 
weak elliptical sources with a minor FWHM larger than 2.0', i.e. 
a full extent along the minor axis larger than ~ 4.0', lost more 
than 15% of their peak flux density and size because of the cor- 
related noise removal of the data reduction. Filaments 1, 3, and 
the widest filament of the T-shaped structure 4 may therefore be 
affected by this spatial filtering and their intrinsic widths may 
be somewhat larger than measured here. On the other hand, the 
measured width of filaments 2, 5, and of the thinnest filament of 
structure 4 likely corresponds to their actual minor size. The to- 
tal mass of these filaments was measured in the full map (Fig.|2]l 
within the contour defined above in the smoothed map at scale 
5. They are listed in Col. 16 of Table [TOl as well as the corre- 
sponding mass per unit length along the major axis in Col. 17. 

With a total flux in the smoothed map at scale 5 of 59 Jy, 
corresponding to a mass of 3 1 M©, these filaments contribute to 
52% of the total 870 fim flux measured in Cha I. Their mean H2 
column density is relatively low, varying from 2 to 5 x 10^' cm"^, 
equivalent to a visual extinction Ay from 2 to 5 mag. When the 
mass included in the smaller scales is added (i.e. in the sum map 
at scale 5), the total mass of these filaments becomes 53 M©, i.e. 
about 88% of the total mass traced with LABOCA in Cha I. 

6.4. Comparison with extinction and 1-0 emission 

6.4.L Correlation between continuum emission and 

The 870 fim dust continuum emission map of Cha I and the 
extinction map derived from 2MASS are overlaid in Fig. [12] 
The overall correspondence is relatively good, especially in field 
Cha-North. In the other fields, the extinction peaks are somewhat 
shifted compared to the 870 fim peaks, but these shifts (< 2') 
are smaller than the angular resolution of the extinction map 
and may not be significant. The only exception is close to Cha- 
MMS 1 (C 1 ) in field Cha-Center, where the extinction map peaks 
about 5' north-east of the 870 yum peak. This offset still remains 
when the 870 fim map is smoothed to 3', showing that it does not 
result from the different angular resolutions of the original data. 
An explanation could be that the extinction map does not trace 
well the H2 column density in this high-extinction region. The 
average intrinsic colors of the small number of background stars 
detected in each resolution element of this part of the map may 
deviate significantly from the assumed average intrinsic colors 
(Sect. 12.2b . Another contributing effect could be the spatial fil- 



tering of the 870 yum data. The extinction map (and the C'^O 1-0 
map, see Sect. I6.4.2l i suggests a structure more extended toward 
the north-east than traced by LABOCA. This missing extended 
structure could bias the peak position of the 3'-smoothed 870 yum 
map toward the south-west. 

The Act H2 column density sensitivity limit of the 870 /urn 
map is 4.2 x 10^' cm~^, which corresponds to Ay ~ 4.5 mag. 
Most of the continuum emission detected at 870 fim is above the 
contour level Ay = 9 mag, and most extended regions with 4.5 < 
Av < 9 mag traced by the extinction map are not detected at 
870 fim. This is due to the spatial filtering related to the sky noise 
removal (see Appendix lA. lb since these low-extinction regions 
have sizes on the order of 10'-20', larger than the field of view 
of LABOCA, especially in fields Cha-Center and Cha-South. On 
the other hand, there are also a few 870 yum clumps detected 
between 4.5 and 9 mag that are not seen in the extinction map, 
most likely because of its poor angular resolution (e.g. Chal- 
C13 and C22 in the southern part and Chal-C20 and C45 in 
the eastern part). A few YSO candidates are even found below 
3 mag (e.g. Chal-S5, S17, and S18 in field Cha- West, Chal-S20 
and S21 in field Cha-North). 

6.4.2. Correlation between 1-0 emission and Ay 

The extinction map derive d from 2MASS is com pared with the 
C'^O 1-0 map obtained by iHaikala etaTI (l2005h with the SEST 
in Fig.[T3h. We smoothed the C'^O 1-0 map to an effective angu- 
lar resolution of 3' and a pixel size of 1.5', to match those of the 
extinction map. The coiTelation between the two tracers looks 
very good above Ay = 6 mag where the C''^0 1-0 integrated 
intensity is proportional to the visual extinction in good approx- 
imation (see Fig. [T3b). The high-extinction end (Ay > 15 mag) 
is however not well constrained. 

At low extinction (Ay < 6 mag), the C'^O 1-0 integrated in- 
tensity drops steeply with decreasing extinction. This sharp de- 
crease is most lik ely due to photodiss ociation. The photodisso- 
ciation models of Visser et al ] (I2009h show such a steep C'^O 
column density variation in the range of extinction they investi- 
gated (slope ~ 2 to 7 in a log-log diagram for Ay < 5 mag, de- 
pending on the density and the UV field). The extension of their 
model to Ay ~ 15 mag (Visser 2010, priv. comm.) shows that for 
a density hh = 300 cm"-', the turnover occurs at Ay ~ 4 - 5 mag 
for a UV-field scaling factor of 1 - 10 and a gas temperature of 
15 K, while it occurs at Ay ~ 2 - 3 mag for = 1000 cm"^ 
and at Ay ~ 7-8 mag for tin = 100 cm"^. The density of 
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Fig. 11. Mass distribution dN/dM (a) and dA^/dlog(M) (b) of 
the 60 starless (or Class 0) sources. The error bars represent the 
Poisson noise (in Va^). The vertical dotted line is the estimated 
90% completeness limit. The thick solid line is the best power- 
law fit performed on the mass bins above the completeness limit. 
The best fit exponent, a and aiog, respectively, is given in the up- 
per right corn er of each pan el. The IMF of single stars corrected 
for binaries (Kroupa 2001, KOI) and the IMF of multiple sys- 
tems (Chabrier 2005, COS) are shown in dashed (red) and dot- 
dashed (blue) lines, respectively. They are both vertically shifted 
to the same number at 5 M©. The dotted (purple) curve is the 
typica l mass spectrum of CO clumps (Blitzlll993t iKramer et alj 
[l998h . 



free particles corresponding to the median extinction in Cha I 
is ~ 350 cm"-' (see Sect. I4.2l i. i.e. hh ~ 590 cm"-'. This den- 
sity, which should characterize the density close to the edges of 
the cloud better than the mean density, should imply a turnover 
be tween Ay ~ 2-3 and 4-5 mag according to the models 
of IVisser et al.L A reason why the turnover in Cha I occurs at a 
higher level (Ay ~ 6 ma j^ may be that the low-density medium 
is inhomogeneous and the UV radiation can penetrate deeper in 
the cloud than predicted by the homogeneous model of Visser et 
al.. 

It follows from the previous paragraphs that the C'**0 1-0 
emission does not trace well the cloud mass at Ay < 6 mag but 
that it is well correlated with the extinction between Ay - 6 and 
~ 1 5 mag and should thus be a reliable tracer of mass in this range 



11 00 00 
a (J2000) 



10 50 00 



Since this visual extinction is about the threshold above which Ry 
in Cha I is thought to have a value higher than the standard one (see 
Appendix IB.3t , we may have to apply a correction to this visual ex- 
tinction that was computed assuming the standard extinction law. With 
Ry = 5.5, the true visual extinction would become A^ ~ 5 mag instead 
of 6, still high compared to the model predictions. 



Fig. 12. Extinction map of Fig. [T] (red contours) overlaid on the 
870 jum continuum emission map of Cha I (black contours). The 
contour levels of the extinction map increase from 3 to 1 8 mag 
by step of 3 mag. The thicker red contour corresponds to Ay = 
6 mag. The contour levels of the 870 //m map are the same as 
in Fig. |2] The dotted line delimits the field mapped at 870 /im. 
The field of view of LABOCA and the angular resolution of the 
extinction map are shown in the upper right corner. 

of extinction. We would then expect C'^O 1-0 to trace about the 
same mass as the extinction map ab ove Ay - 6 mag (220 Mq, 
see Sect. gJli. iHaikala et afl (l2()05h derived a mass of 230 Mq 
from the C'^O 1-0 map shown in Fig. [T4k . This map d oes not 
contain field Cha- West. This field was covered by Mizun oet al.l 
(1199 9) and its total C'^O 1-0 integrated intensity represents an 
additional contribution of 15%. The mass traced by C'^O 1-0 
in Cha I is thus about 260 Mq. Given the uncertainties on the 
C'^O abundance, the excitation temperature, and the Ay-to-mass 
conversion, this is in very good agreement with the mass derived 
from the extinction map above Ay - 6 mag. 

6.4.3. Comparison of continuum emission witli C"*0 1-0 
emission 

The 870 fim dust continuum emission is much less correlated 
with the C'^O 1-0 emission than the 2MASS extinction, as 
shown in Fig. [T4h . We smoothed the 870 fim continuum map 
to an effective angular resolution of 45" and a pixel size of 30" 
to match the C'^O map resolution and pixel siz^El The ratio 
of C'^O 1-0 integrated intensity to 870 fj.m peak flux density is 
shown in Fig.fT4b in color scale. It varies by more than one order 
of magnitude, being the lowest at the 870 yum peaks. If we as- 
sume that the dust temperature and the excitation temperature of 



" But beware that the C"*0 map is highly undersampled since it was 
observed with a step of 60" . 
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Fig. 14. a Integrated intensity map of C"*0 1-0 emission detected with the SEST bv lHaikala et alJ (l2005l) in T* scale. The 870 yum 
continuum map smoothed to 45" is overlaid as contours. The contours are -/ (dashed), /, 2/, 4/, 8/, 16/, and 32/, with / = 4cr 
and (T the rms noise level (23 mJy/45"-beam). The dotted line delimits the field mapped at 870 jim. The filled circle shows the SEST 
angular resolution, and the square the sampling of the C'**0 map. b Ratio of C'**0 1-0 integrated intensity to 870 jjm continuum 
flux density with a 870 jjm clip at Act. The contours are the same as in a. c Correlation between the C"*0 1-0 integrated intensity 
and the 870 |im continuum emission. Only data points with a C'**0 1-0 integrated intensity higher than 0.2 K km s"' (> 3cr for 90% 
of the map) and a 870 /im flux density higher than 69 mJy/45"-beam (3cr) are plotted. The red line shows the result of a linear fit 
over the full range. The blue dashed line has a slope of 1 in log-log space and its vertical position is arbitrary. 



C'**0 1-0 are the same, the g as-to-dust mass ratio is uniform, and 
the line is optically thin (see iHaikala et al.ll2005l) . then the ratio 
plotted in Fig. [T4b should roughly trace the CO abundance. Its 
variations point to significant CO depletion toward the center of 
the 870 fim sourceo. Figure [T4b displays the relation between 
the C'^O 1-0 integrated intensity and the 870 ixm flux density 
for all pixels above 3cr. Again, this plot shows the poor corre- 
lation between both tracers. In addition, a linear fit in log-log 
space yields a slope of ~ 0.2, far from what would be expected 
if C'**0 1-0 were tracing the 870 ;um continuum emission well 
(see the red versus blue lines in Fig.[T4b). 

However, these figures are a priori biased because of the spa- 
tial filtering of the 870 yum data that does not affect the C''^0 1-0 
data. The mass traced by the 870 /im map represents only ~ 30% 
of the mass derived from C'^'O 1-0 (see Sect. l4.2l i. Since the esti- 
mate of the latter does not take into account depletion in regions 
of high column density, the discrepancy between both tracers 
must be even worse than a factor of 3. A proper correction would 
be to filter the C'**0 1-0 data in the same way as the 870 jum data 
are filtered. This is not an easy task and, because of the under- 



The temperature dependences of the C'*0 and H2 column densities 
when derived from the C'^0 1-0 integrated intensity and the 870 yum 
flux density are such that their ratio decreases by a factor 5 when the 
temperature drops from 12 to 6 K. Therefore, if the temperature is not 
uniform but rather decreases toward the center of each core, the drop of 
CO abundance toward the center of each core would be even more se- 
vere than derived above with the assumption of a uniform temperature. 



sampling of the C'^O 1-0 data, may not be fully reliable either. 
At this stage of the analysis, we cannot exclude that a significant 
part of the spatial variations seen in Fig.[T4b are due to the filter- 
ing of the continuum data rather than solely to variations of the 
C'**0 abundance. 



7. Discussion 

7.1. Evolutionary state of Cha-MMS1 and Ciia-MMS2 

Cha-MMSl and Cha-MMS2 were the only tw o dense cores 
known in Cha 1 before our LABOCA survey. iReipurth e t al.l 
(Il996h discovered them at 1.3 mm with the SEST and sug- 
gested that both could be Class protostars based on their tenta- 
tive association with known outflows. In Cha-North, the peak of 
the 870 jum emission detected with LABOCA (s ource S3) cor- 
respon ds to the source Cha-MMS2 detected by IReipurth et al.l 
(119961) within 7", which is certainly within the errors of the 
SEST measurements that were performed with a single-channel 
bolometer. Our source S3 is clearly associated with the T Tauri 
star WW Cha (coincidence within 1.2", see Table |2|. As a re- 
sult, Cha-MMS2 is associated with WW Cha and i s certainly 
not a Class protostar (see also iLommen et al ] l2009l who came 
to the same conclusion). 

Belloche et al. (2006) discovered a faint Spitzer source at 
24 and 70 pm located 5-7" to the east of Cha-MMSl (source 
A^sp = 1 in Table |4|. They also reported a Spitzer non-detection 
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Fig. 13. a Integrated intensity map of C "^ 1-0 emission de- 
tected with the SEST by Hai kala et all (l2005h in T* scale, 
smoothed to an effective angular resolution of 3' (filled circle). 
The extinction map of Fig.[T]is overlaid as contours that increase 
from 3 to 18 mag by step of 3 mag. b Correlation between the 
C'^O 1-0 integrated intensity and the visual extinction. Only 
data points with a C'^O 1-0 integrated intensity higher than 0.1 
K km s~' (> 3cr for 90% of the map) and Ay > 2 mag are plot- 
ted. The dotted line marks the extinction Ay - 6 mag above 
which a linear fit with a fixed slope of 1 (in log-log space) was 
performed. The red line shows the result of this fit. Note that 
all visual extinction values where derived assuming the standard 
extinction law and have to be multiplied by 0.85 for Ay above 
~ 6 mag (see Appendix lB.3l i. 



at 8 j-im (see also Fig. 8 of iLuhmanI l2008l) . at a 3cr level ~ 9 
and 7 times lower than the peak fluxes of the Taurus Class 
protostars IRAM 04191 + 1522 and LI 52 IF, respectively. They 
did not find any evidence for an outflow associated with Cha- 
MMSl based on CO 3-2 observations done with APEX, and 
derived a deuterium fractionation of 1 1% for NaH^, a value typ- 
ical for evolved prestellar cores and young protostars. Based on 
these two results, they concluded that Cha-MMS 1 could be at the 
stage of the first hydrostatic core, i.e. just before the dissociation 
of H2 and the second collapse. Our LABOCA observations re- 
fine the position of Cha-MMS 1 . The 870 jim emission peaks at a 
position 1.3" away from the Spitzer source, which is well within 
the positional accuracy of our measurement. This strengthens 
even furflier the association of Cha-MMS 1 with the faint Spitzer 

source. 

Recently, i Dunham et alJ (l2008h found a very tight correla- 
tion between the Spitzer 70 pm flux and the internal luminos- 
ity of the low-luminosity YSOs of the c2d legacy project. They 
define the internal luminosity as the sum of the luminosities 
of the central protostar and its circumstellar disk, if present. 
With a 70 pm flux of ~ 200 mJy (Belloche et al. 2006), this 
correlation implies that Cha-MMS 1 has a very low internal lu- 
minosity of ~ 0.015 Lq. Luminosities ranging between 10""* 
and 0.1 Lq have been theoretically predicted at the stage of 
the first core, depending on the collapse model and the amount 
of rotation (Saigo & Tomisaka 2006; Omukai 2007). The faint 
luminosity of Cha-MMS 1 looks therefore consistent with the 
luminosity of a first core. It is about 2 and 3 times fainter 



than the internal luminosities of the young Cla ss protostars 
L152 1F and IRAM 04191-H1522, respectively (Dunham et alj 
I2OO8I). Sin ce IRAM 04 1 9 1 + 1 522 drives a prominent bipolar out- 
flow (Andre et al. 1999) and LI 52 IF may also drive a very com- 
pact one (Shinnag a et ani2009l) . the absence of evidence for an 
outflow driven by Cha-MMS 1 is an additional argument suggest- 
ing that it is in a different stage than these two young protostars, 
i.e. that it may be a t the stage o f the fi rst core. 

We note that iChen et all (1201 Oh proposed L1448 IRS2E 
as a candidate first hydrostatic core based on its low bolo- 
metric luminosity (Lboi < 0.1 L©, with no Spitzer detection 
at 70 pm and below) and their detection of a compact out- 
flow. Magnetohydrodynamical simulations do predict that a first 
core may drive an outflow but they predict that it should be 
slow (a couple of km s ' at most, see, e.g.. 'Tom isaka! I2OO2I : 
Machidaet al. 2008|). The outflow of L1448 IRS2E has veloc- 
ities about one order of magnitude higher than the predictions 
of Machida et al.. (2008), more consistent with velocities pre- 
dicted for the second outflow produced once the protostar (sec- 
ond core) is formed. It is therefore somewhat unclear whether 
L1448 IRS2E can really be at the stage of the first core. 

On the other hand, the absence of evidence for an outflow 
powered by Cha-MMS 1 does not contradict these MHD simu- 
lations since a compact, slow outflow is a prior i not ruled out 
in this source. We thus confirm the suggestion of iBelloche et al.l 
(2006) that Cha-MMS 1 is a candidate first hydrostatic core. 

7.2. Nature of the LABOCA starless sources 

In the past decade, the mass distribution of starless dense cores 
in many nearby molecular clouds as traced mainly by (sub)mm 
continuum emission (see references in Sect. 16.1.71 ) was found 
to have the same power-law behaviour as the stellar IMF to- 
ward the high-mass end, while the mass distribution of lower- 
density clumps traced with CO has a much flatter shape. This 
was put forward as a strong argument suggesting that the IMF 
in these low-mass-star forming clouds is essentially set at the 
early stage of cloud fragmentation already. The consistency of 
the high-mass end of the mass distribution of the Cha I starless 
sources with the stellar IMF suggests that they may be prestel- 
la^in nature, i.e. form stars in the future. 

However, the distribution deviates significantly from the stel- 
lar IMF at th e low-mass en d. This CMD is compared to the 
st eUar IMF of iKroupal (12001) and the IMF of multiple systems 
of lChabrieH (l2005h in Fig. [TT] Since LABOCA does not probe 
the scales at which members of a (future) multiple system can 
be seen individually, only the comparison to the IMF of mul- 
tiple systems should make sense. While the Cha I mass distri- 
bution looks consistent with the IMF of multiple systems above 
~ 0.3 Mq, it deviates significantly below this threshold, the Cha I 
starless sources being more numerous at the low-mass end than 
expected from the IMF. Since the completeness limit of the Cha I 
survey is ~ 0.2 M©, the discrepancy is a priori even worse at 
very low mas s. Such a de viation has already been noticed in 
Ophiuchus bv I Andre et al.l (12007) who concluded that a fraction 
of the lowest-mass starless sources may not be prestellar and that 
the link between the CMD and the stellar IMF is l ess robust at 
the low-mass end than at the high-mass end (see also lAndre et al.l 
2009). Our LABOCA survey of Cha I supports this conclusion. 
Follow-up heterodyne observations with MOPRA and APEX to 
measure the fraction of gravitationally bound sources are under- 



'-^ A prestella r core is usually defined as a starless core that is gravita- 
tionally bound l lAndre et al.ll2000l:IWard-ThomDson et al.ll2007h . 
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way and should help better understanding the low-mass end of 
the mass distribution in Cha I. 

Only four sources - among them Chal-Cl, i.e. Cha-MMSl, 
that is not starless (see Sect. 17. Il l - have a mass larger than the 
critical Bonnor-Ebert value, and seven additional ones may also 
belong to this group if their mass is underestimated by a fac- 
tor of 2 (see Sect. I6.I.6I 1. As a result, at most 17% of the star- 
less sources appear to be gravitationally unstable and thus likely 
prestellar, provided their turbulent and magnetic supports are 
small compared to their thermal support or will vanish if they 
are not. This fraction^is s ignific antly smaller than in Ophiuchus 
and Aquila. lAndre et al] (i20i0al) estimate that more than 60% of 
the 541 starless sources detected with Herschel in Aquila, with a 
completeness limit similar to ours (see Sect. 16.1. St . are candidate 
prestellar cores (with MfiE/Mt ot < 2). Sim ilarly, about 70% of 
the starless sources detected bv lMotte et al. ( 1998) in Ophiuchus 
with a completeness limit somewhat better than ours have 
MBE(R)/Mtot < 2 (see definition of Mbe(R) in Sect. |6X6] R 
In Taurus, more t han 5 0% of the starless cores identified by 
iKauffmann et al.l (12008) are gravitationally unstable, but this 
sample is biased and inhomogeneous so a direct comparison to 
Cha I is more uncertain. 

The very low fraction of candidate prestellar cores among 
the starless cores detected in Cha I suggests that we may be wit- 
nessing the end of star formation in this cloud or, on the contrary, 
that these starless cores are in a less advanced stage of evolution 
compared to those in Aquila, Ophiuchus, and Taurus. The latter 
suggestion is less likely because we already assessed the gravi- 
tational state of the cores assuming that they have no turbulent 
and magnetic support. A reduction of these support mechanisms 
via turbulence dissipation or ambipolar diffusion could in prin- 
ciple lead to gravitational instability and eventually gravitational 
collapse, but only if the cores are above the (purely thermal) crit- 
ical Bonnor-Ebert mass limit already. A reduction of the thermal 
support via cooling could in principle also lead to a gravitational 
instability. In Sect. 16.1.6! we estimated that even if the current 
temperature of the cores is 7 K instead of 12 K, this would not 
change the current fraction of unstable cores significantly. If the 
cores are currently at 12 K even in their center (i.e. the masses 
we measure are correct) and are still cooling down, then the frac- 
tion of unstable cores will not change much either, because the 
Bonnor-Ebert critical mass will decrease by less than a factor of 
2 while most cores are at least a factor of 3 below the current 
critical mass limit. It is thus unlikely that most starless cores that 
appear to be currently stable will become unstable in the future, 
unless they experience a strong external perturbation. 

7.3. Star formation efficiency 

7.3.1 . Local star formation efficiency 

Apart from the deviation at low mass, the mass scale of the Cha I 
CMD agrees quite well with that of the IMF of multiple systems 
in Fig. [TT] Taken at face value, this would suggest that the star 
formation efficiency at the core level in Cha I will be as high as 
~ 100% for all prestellar sources t hat populate t he high mass end 
above ~ 0.3 Mq. For comparison. I Andre et al.l (12010 a') found an 
effici ency fcoie = Mi,/Mcore ~ 20^0% in Aquila and Alves et al. 
(2007) obtained Ecok - 30± 10% for the Pipe nebula. Enoch et alj 
t2008) could only set a lower limit of 25% for their combined 



The percentage quoted here assumes R = FWHM and a distance 
of 125 pc. 



sample of Perseus, Serpens, and Ophiuchus sources, because of 
their poor completeness limit (0.8 Mq). 

Luhman (2007) computes the specific IMF of the stellar pop- 
ulation of Cha I and find a peak at 0.1-0.15 Mq. If we take this 
turnover mass as reference, then the (apparent) turnover of the 
Cha I CMD at the completeness limit still suggests a local star 
formation efficiency ecore > 50-70%. Such a very high efficiency 
for Cha I sounds a priori unlikely. One explanation could be that 
the contamination by gravitationally unbound, starless sources 
(see Sect. I7.2l i may extend further to higher ma sses and mask 
the maximum of the true prestellar CMD. Indeed. lKonvves et al.l 
(2010) found that the peak of the Aquila CMD shifted from 0.6 
to 0.9 Mq when they removed from thek sample the sources that 
they estimated not to be prestellar (~ 35%). On the other hand, 
the majority of the Pipe nebula sources appear to be gravitation- 
ally unbound (Lada et al. 2008) and th eir CMD still closely re- 
sembles the IMF in shape dAlves et al.' '20b7). Therefore it is un- 
clear by how much the Cha I CMD could be affected by contam- 
ination from unbound, starless sources. The on-going follow-up 
heterodyne observations should help answer this question. 

7.3.2. Global star formation efficiency 

At the scale of the cloud, iLuhmanI (l2007h estimate a current 
star formation efficiency, defined as SEE - jf^^-^, of about 
10%, based on the cloud mass derived from CO (~ 1000 Mq, 
see Sect. I4.2l i and the mass ~ 100 Mq of the 215 young stars 
and brown dwarfs known to belong to Cha I at that time. This 
star formation efficiency is a factor of 1.5 to 3 larger than the 
SFEs of the five nearby molecular clouds targeted by the c2d 
Spitzer legacy project (between 3.0 and 6.3%, see Table 4 of 
Evans et al. 2009). It is even more than one order of magnitude 
higher than the SEE of the Taurus molecular cloud (0.6 see 
Goldsmit h et al. 200 8). Assuming a constant rate of star forma- 
tion, Ev ans et al.l (12009!) suggest that the c2d clouds could reach 
a final SEE of 15 to 30% given the estimated typical lifetime of a 
molecular cloud. The SEE of Cha I is close to this limit, this be- 
ing even more true if the rate of star formation in the c2d clouds 
decreases with time like in Cha I (see Sect. l7.4l i. which would re- 
duce the expected final SEE. We note that Bon nell et al.l (1201 Ol) 
find an overall efficiency of star formation of ~ 15% in their 
numerical simulation of a turbulent, molecular cloud that un- 
dergoes star formation. It is thus very likely that Cha I is much 
closer to the end of the star formation process than the other 
clouds mentioned above. 

7.4. Lifetimes 

The median age of the YSO population of Cha I is ~ 2 Myr 
(Luhman 2008|3- With the census of starless cores and YSOs 
presented in Table (TT] we can statistically estimate the lifetimes 
of each category if we assume that all prestellar cores evolve 
into Class 0/Class I protostars and then into Class Il/Class III 
YSOs, and that star formation has proceeded at a constant rate 
over time. Five to 15 YSOs out of 204 are protostars, implying a 
typical lifetime of 0.5 - 1.5 x 10^ yr for the protostars in Cha I. 
Only one object is a Class protostar (or first core), hence a 



This SFE was estimated with a cloud mass that includes very low 
column-density material. If we take the mass of the regions with Ay > 
3 mag (see Table ll2! (. then the SFE would be ~ 4 times larger. 

If the further distance is confirmed for Cha I (196 pc, see 
Appe ndix IB. il l, the median age would be reduced to ~ 1 Myr dKnudd 

Mm- 
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Table 11. Census of starless cores and YSOs in Cha I. 



Class Number 

Starless cores 59° 

Class (or first core) 1* 

I 4' 
Flat 10'' 

II 94' 

III 95^^ 



Notes. Th i s wor k. <''' Cha-MMSl iBelloche et al]l2006h . <'^' From 
iLuhman et aO (l2008h . over a large (but not complete) part of Cha I. 

very short lifetime of 10^ yr for the Class phase. Given that a 
fraction of the starless cores may not be prestellar in nature (see 
Sect. l7.2l i. we can only derive an upper limit of 6 x 10^ yr to the 
duration of the prestellar phase in Cha I. 

The median peak density of the Cha I starless cores is 
3.5 X 10^ cm"^ and their median mean density within 50" is 
7 X 10"* cm"^ (Figs. [8^ and f). This median mean density corre- 
sponds to a prestellar lifeti me of ~ 4x 10^ yr acco rding to the em- 
pirical relation of iWard-Thompson et al.l (l2007h . If the fraction 
of non-prestellar sources among the Cha I starless sources is on 
the order of 30% or less, then the agreement between this prestel- 
lar lifetime and the lifetime derived from the statistical counts in 
the previous paragraph suggests that the assumption of a star for- 
mation rate constant over time until now is valid in Cha 1. It is 
also consiste nt with the prestella r lifetime of 4.5 + 0.8 x 10^ yr 
estimated by lEnoch et alJ (l2008h for their combined sample of 
Perseus, Serpens, and Ophiuchus (see also |j0rgensen et alj|2007l 
for Perseus alone). The very short lifetime of 10"* yr derived here 
for the Class phas e is similar to the one found in Ophiuchus 
(1 - 3 X 10"* yr, see I Andre et al.ll2000i). but about orie orde r of 
magnitude shorter than in Perseus (e.g. lHatchell et al1l2007h . 

On the other hand, if, e.g., 90% of the Cha I starless sources 
turn out not to be prestellar in nature, then the unrealistically 
low statistical lifetime that would be inferred for the prestellar 
sources (6 x 10"* yr) would rather point to a decrease of the star 
formation rate with time and confirm the trend pointed out by 
iLuhmanI (120071) . This could also explain why so few Class 
sources (only one) are found in Cha 1, if the Class lifetime 
derived in Perseus rather than in Ophiuchus applies to Cha I. In 
any case, an increase of the star formation rate until the present 
time (.Palla & Stabler 2000. see Sect.[B is ruled out for Cha I. 

7.5. Filaments 

The Herschel satellite has recently revealed an impressive 
network of parsec-scale filamentary struc tures in the nearby 
molecular clouds Aquila and Polaris (e.g. lAndre et"aD 1201 Oat 
en'shchikov et al. The few filaments detected in Cha I 

with LABOCA (see Sect. [Oil and in C"*0 1-0 (Haika la et all 
I2005h are strong evidence for the existence of such a net- 
work in this cloud, of which we only see the tip. These fila- 
ments have a deconvolved width of 0.16 + 0.05 pc on average 
(33000 AU), measured at a level of 10 mJy/21.2"-beam. This 
width may be somewhat underestimated because of the spatial 
filtering due to the correlated noise removal (see Sect. 16.3b . It 
is similar to the typical width o f the Perseus filaments as mea- 
sured by iHatchell etaLl (l2005h with SCUBA {FWHM ~ 1', 
i.e. 15000 AU at 250 pc). However, the Cha I filaments are 
less massive, with a mass per unit length ranging between 8 
and 33 Mq pc ' (see Table [TOb while i t ranges between 47 and 
1 15 Mq pc ' for the Perseus filaments dHatcheU et al.ll2005h . 



The mass per unit length of filaments is an essential parame- 
ter to evaluate their stability. In Aquila, most (> 60%) candidate 
prestellar cores were found in (thermally) su percritical filaments 
defined by Mn„ > M™™f ^ la^/G (see lAndre et al.l 1201 Oai 
and references therein). Only one filament in Cha I (Nm = 1) has 
a mass per unit length larger than the critical value (19 pc"' 
at 12 K). It contains the single Class protostar (Chal-Cl) and 
two of the three candidate prestellar cores (Chal-C2 and C3), 
while only one of the four other (subcritical) filaments con- 
tains a candidate prestellar core (Chal-C4 in filament 5). Even 
if these numbers are much too small to derive any robust con- 
clusion, they seem to be consistent with the results obtained by 
lAndre et all (|20 10a) for Aquila. 

However, the total velocity dispersion as probed by 
C'**0 1-0 in Cha I is significantly larger than the purely ther- 
mal dispersion. The C'**0 1-0 linewidth (FWH M) is typically 
on th e order of 0.5-0.8 km s ' (see Table 1 of iHaikala et al.l 
2005"), yielding a non-thermal velocity dispersion cr nt — 0.20- 
0.38 km s"' and an eff'ective sound speed Oeff = 0.29- 
0.43 km s"' at 12 K. We can compute the mean free-particle den- 
sity nmean of the filaments with the equation nmean - ^t'^ x 
with A^inean and a„iin given in Table [10] and jju^ and yu de- 
fined in Appendix IB. 5 1 These mean densities are in the range 
3.7-10 X 10^ cm"^ Since the critical density of the C"*0 1-0 
line at 10 K is about 3 x 10^ cm"^, this transition traces the 
bulk of the filaments very well and its linewidth should reli- 
ably trace their total velocity dispersion, although some con- 
tamination from the ambient cloud cannot be completely ex- 
cluded. If we assume that C''^0 1-0 traces only the filaments, 
then the critical mass p er unit length become s the virial mass 
per unit length 2 a^^ /G (iFiege & Pudritz|[2000l) and ranges from 
39 to 86 Mq pc The Cha I filaments have a mass per unit 
length lower than this critical value by a factor of 2 at least 
(the velocity dispersion is the highest in filament Nfi\ - 1). 
Unles s they are compresse d by a significant toroidal magnetic 
field (iFiege & Pudritz"2000) or the C'**0 1-0 Hnewidth overes- 
timates their total velocity dispersion, they do not seem to be 
supercritical and will most likely neither collapse nor fragment 
in the future. This would be consistent with our conclusion that 
the process of star formation is nearly over in Cha I. 

7.6. The end of star formation in Cha 1 ? 

The low fraction of candidate prestellar cores among the popula- 
tion of starless cores detected with LABOCA, the small number 
of Class protostars, the high global star formation efficiency, 
the decrease of the star formation rate with time, and the low 
mass per unit length of the detected filaments all suggest that we 
may be witnessing the e nd of the s t ar fo rmation process in Cha I. 
While the conclusion of iLuhmaiil (1200 7) that the star formation 
rate in Cha I has decUned fo r several Myr already questions the 
scenario o f iPalla & Stable^ ("2000) in which star formation is a 
threshold phenomenon and accelerates as a result of cloud con- 
traction (see Sect. [1]), o ur findings raise new questions about the 
end of star formation. Palla & Stahleil (l2000l) argue that the trun- 
cation period, i.e. the period over which the star formation rate 
drops, must be very brief (1 Myr at most) and that the falloff 
should be related to the dispersal of the parent cloud gas. In low- 
mass star forming regions, they consider that this dispersal can 
only result from the impact of molecular outflows. Even if the 
total cuiTent mass efflux due to molecular outflows in the partic- 
ular case of Taurus is three orders of magnitude lower than the 
one needed to disperse the cloud in such a short time, they argue 
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Table 12. Mass, area, and average surface density of a sample of nearby molecular clouds, measured over the regions with Ay > 
3 mag, 3 < Av < 6 mag, and Ai/ > 6 mag. 









Ay > 3 mag 


3 < Av < 6 mag 


A 


V > 6 


mag 








Cloud 


Distance 


Mass 


Area E" 


Mass 


Area 




Mass 


Area 


E" 


Rm' 


Ra" 


R,." 




(pc) 


(Mo) 


(deg^) (Mo pc-^) 


(Mo) 


(deg-) 


(Mo pc--) 


(Mo) 


(deg^) 


(Mo pc--) 








(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Chamaeleon I 


150 


720 


1.2 88 


500 


0.9 


81 


220 


0.30 


109 


0.45 


0.33 


1.4 


Ophiuchus 


125 


3340 


10.1 69 


2350 


8.3 


60 


990 


1.87 


111 


0.42 


0.23 


1.9 


Taurus 


140 


6120 


15.1 68 


4320 


12.5 


58 


1800 


2.59 


116 


0.42 


0.21 


2.0 


Perseus 


250 


6140 


4.8 68 


4510 


3.9 


61 


1620 


0.87 


98 


0.36 


0.22 


1.6 


Serpens 


260 


7980 


5.5 70 


6300 


4.5 


68 


1680 


1.00 


82 


0.27 


0.22 


1.2 


Orion 


450 


63320 


15.0 68 


44180 


11.8 


61 


19140 


3.23 


96 


0.43 


0.27 


1.6 



Notes. The masses and areas for all clouds but Chamaeleon I were provided by S. Sadavoy (priv. comm.), based on lSadavov et alj ( l201(]h . The 
masses for Ay > 6 mag, initially calibrated for Ry = 3. 1, are here reca l ibrated for Ry = 5.5 (see Appendix lB.3t . The masses for Ay > 3 mag are 
corrected accordingly. The distances are the same as in lSadavov et alj d2010h . Average surface density, i.e. mass-to-area ratio. **' Mass, area, 
and surface-density ratios between the regions with Ay > 6 mag and 3 < Av < 6 mag. 



that the (potentially) accelerating nature of star formation should 
significantly increase the mass efflux in the future and that, if star 
formation is really a threshold phenomenon, its falloff can occur 
before the actual cloud dispersal. 

To look for signs of cloud dispersal in Cha I, we compare its 
mass, area, and surface density to those of a sample of nearby 
molecular clouds in Table [121 We list the surface densities for 
the low- (3 < Ay < 6 mag) and high-extinction (Ay > 6 mag) 
regions in Cols. 8 and 11, respectively. The limit Ay = 6 mag ap- 
proximately coiTesponds to the extinction threshold for star for- 
mation mentioned in Sects. [nand l6. 1.11 If all clouds have similar 
depths, then the surface density scales with the average density 
and could be a good indication for cloud dispersal. On the one 
hand, Cha I has a surface density above Ay = 6 mag very sim- 
ilar to the ones of the other clouds, suggesting that the average 
density of its high-extinction region is the same as in the other 
clouds, provided its depth is similar to the ones of these clouds. 
Only Serpens has a significantly lower surface density compared 
to the rest of the sample. On the other hand, the surface density 
of Cha I (and Serpens to a lesser extent too) for 3 < Ay < 6 mag 
is significantly higher than in the other clouds. As a result, it 
stands out, along with Serpens, with a surface-density ratio be- 
tween the regions with Ay > 6 mag and 3 < Ay < 6 mag some- 
what lower than in the other clouds. Since this ratio does not de- 
pend on the depth of the cloud, it could be a more reliable tracer 
of cloud dispersal and could indicate that Cha I and Serpens are 
in the process of being dispersed. However, these two clouds do 
not seem to differ significantly from other nearby star-forming 
clou ds in terms of probabilit y density functions of column den- 
sity dKainulainen et al.ll2009l) . Therefore the signs for cloud dis- 
persal in Cha I based on the extinction map are very tenuous. 
Further clues to explain why Cha I could be at the end of the star 
formation process should be searched for. Finally, we note that 
we cannot exclude that the star formation activity in Cha I could 
be revived in the future if the cloud would undergo a significant 
external compression 



8. Conclusions 

We performed a deep, unbiased, 870 jum dust continuum sur- 
vey for starless and protostellar cores in Chamaeleon I with the 
bolometer array LABOCA at APEX. The resulting 1 .6 deg^ map 
was compared with maps of dust extinction and C'^0 1-0 emis- 
sion. It was decomposed into multiple scales to disentangle the 
emission of compact or moderately extended sources from large- 



scale structures. The analysis was performed by carefully taking 
into account the spatial filtering properties of the data reduction 
process. Our main results and conclusions are the following: 

1 . The standard reduction process of LABOCA data with BoA 
filters the emission more extended than ~ 2.5-4' (FWHM), 
depending on its flux density. Elongated sources are better 
recovered than circular ones along their major axis, up to 6- 
9'. The iterative scheme requires about 21 iterations to con- 
verge. 

2. The mass detected with LABOCA (61 Mq) represents only 
6-8% of the cloud mass traced by the dust extinction and 
CO, and about 30% of the mass traced by C"*0. Most (88%) 
of this 870 fim mass is located in five filaments with typical 
length and width of 0.6 and 0.16 pc, respectively. 

3. The C"^0 1-0 emission is a reliable tracer of mass in the 
visual extinction range 6 to 15 mag. C'^O appears to be pho- 
todissociated below 6 mag. 

4.81 sources were extracted from the 870 jum map, 59 of which 
are starless based on the absence of Spitzer association, one 
is a candidate first hydrostatic core (Cha-MMSl), and 21 are 
associated with more evolved young stellar objects already 
known in the SIMBAD database. 

5. Like in other nearby molecular clouds, the starless cores in 
Cha I are only found above an extinction threshold Ay = 
5 mag and their distribution of extinction peaks at Ay ~ 
8 mag (7 mag with Ry = 5.5). Alternatively there may be 
no strict threshold but just a very low probability of finding 
a core at low extinction. 

6. The 90% completeness limit of our 870 jum survey is 
0.22 M0 for starless cores. Only 50% of the detected star- 
less cores are above this limit, suggesting that we may miss 
a significant fraction of the existing starless cores. 

7. The starless cores in Cha I are less dense than those in 
Perseus, Serpens, Ophiuchus, and Taurus by a factor of a few 
on average. This may result from the better sensitivity of the 
Cha I survey and/or the multi-scale analysis of the Cha I data. 

8. The CO abundance seems to decrease toward the peak of the 
embedded starless cores, as expected for bona-fide prestellar 
cores, but the spatial filtering of the 870 //m data prevents a 
firm conclusion. 

9. On the one hand, the high-mass end of the core mass dis- 
tribution in Cha I is consistent with the stellar initial mass 
function, suggesting that a fraction of the starless cores may 
be prestellar. On the other hand, the low-mass end is over- 
populated, and at most 11% of the starless cores have a mass 
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higher than the critical Bonnor-Ebert mass. Both results sug- 
gest that a large fraction of the starless cores may not be 
prestellar. 

10. The local star formation efficiency at the core level in Cha I 
(> 50%) seems to be larger than in other molecular clouds, 
but this conclusion is uncertain due to possible contamina- 
tion by unbound, starless sources. 

1 1 . The global star formation efficiency at the scale of the cloud 
is 10%. It is significantly larger than in other clouds and sug- 
gests that Cha I is much closer to the end of the star forma- 
tion process. 

12. The star formation rate in Cha I can have been constant over 
time only if > 70% of the detected starless cores turn out to 
be prestellar This is very unlikely, which suggests that the 
star formation rate has decreased over time, as also indicated 
by the very small number of Class protostars. An increase 
of the star formation rate to the present time in Cha I is ruled 
out. 

13. Apart from maybe one, the filaments in which most starless 
cores are embedded do not seem to be supercritical and will 
most likely neither collapse nor fragment in the future. 

The low fraction of candidate prestellar cores among the 
population of starless cores detected with LABOCA, the small 
number of Class protostars, the high global star formation effi- 
ciency, the decrease of the star formation rate with time, and the 
low mass per unit length of the detected filaments all suggest that 
we may be witnessing the end of the star formation process in 
Cha I. We only find a tenuous evidence for cloud dispersal based 
on the distribution of extinction. Further clues to explain why 
Cha I could be at the end of the star formation process should be 
searched for. 
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Appendix A: LABOCA data reduction: spatial 
filtering and convergence 

A.1. Spatial filtering due to the correlated noise removal 

Since LABOCA is a total-power instrument used without wob- 
bler, the standard data reduction relies on the correlated signal 
of all pixels and of smaller groups of pixels to estimate and re- 
move the atmospheric and instrumental noise. This method acts 
as a high-pass filter and makes the detection of very extended 
structures extremely difficult. The strongest limiting factor is the 
size of the groups of pixels co nnected to the same read-out ca- 
ble (about 1.5' by 5', see ISchuUer et alj|2009ft . To analyse the 
dust continuum emission map in detail, it is therefore essential 
to understand the properties of this spatial filtering. 

We performed several Monte Carlo simulations to character- 
ize the filtering. First, we inserted artificial, circular, Gaussian 
sources in the raw time signals with the BoA method addSource 
after calibration. The sources were inserted at random positions 
in signal-free regions of the map shown in Fig. |2] The peak flux 
densities and FWHM sizes of the sources were randomly as- 
signed in the range 50 to 500 mJy/beam and 19.2" to 8', respec- 
tively. The only constraint was that the artificial sources should 
not overlap and should not be too close to the map edges either. 
The modified raw data were then reduced following exactly the 
same procedure as described in Sect. 13.11 except that the final 
map was not smoothed and kept its original angular resolution 
of 19.2" (HPBW). Each artificial source was fitted with an ellip- 
tical Gaussian using the task GAUSS_2D in GREG to measure 
its output peak flux density and sizes after data reduction. We 
also measured the total output flux on a circular aperture of ra- 
dius equal to FWHM. 

The whole procedure was repeated twice, with two indepen- 
dent samples of 30 artificial sources. Figures lATh -e present the 
results. As a control check, we show in Figs. lA.ll F-i the results 
of the fits performed on the artificial sources directly inserted 
in the final unsmoothed continuum emission map of Cha I. The 



Table A.l. Filtering of artificial Gaussian sources due to the cor- 
related noise removal depending on the input size along the mi- 
nor (major) axis (FWHM). 



Source sample 




Size where 




(1) 


/;//; = 85% 

(2) 


sys' = 85% 
(3) 


f;^,/fL = 85% 
(4) 


Circular sources: 








strong sources 

{> 150 mJy/beam) 


280" 


250" 


170" 


weak sources 

{< 150 mJy/beam) 


170" 


140" 


90" 


all sources 


240" 


200" 


120" 


Elliptical sources with aspect ratio 2.5: 

strong sources 220" (550") 220" (550") 

(> 150 mJy/beam) 


130" (330") 


weak sources 

(< 150 mJy/beam) 


130" (330") 


120" (300") 


70" (180") 


all sources 


180" (450") 


160" (400") 


100" (250") 



Notes. f° and are the output and input peak flux densities, s^ and s' 
the output and input sizes (FWHM), and and fl^^^ the total output 
and input fluxes over an aperture of radius equal to the input FWHM. 
The sizes in parentheses correspond to the major axis of the elliptical 
sources. 



effects of the spatial filtering due to the data reduction are ob- 
vious when comparing Figs. lA.lT ^-e to Figs. lA.lt -i. The noise 
produces a small dispersion in the output source properties, es- 
pecially for weak sources (crosses) but apart from that, all curves 
are basically flat in the control sample (Figs. lA.lT -i). On the con- 
trary, after data reduction, the peak flux densities of the sources 
are reduced by ~ 15% to ~ 70% for input sizes varying from 
~ 240" to ~ 440". Similarly, the sizes of the sources are re- 
duced by ~ 15% to ~ 50% for input sizes from ~ 200" to 
~ 440". The total fluxes suffer the most from the spatial fil- 
tering: they are reduced by ~ 15% to ~ 90% for input sizes 
from ~ 120" to ~ 440". A careful inspection of Figs. lA.l^ - 
j shows that the sources with a peak flux density weaker than 
~ 150 mJy/beam (crosses) are systematically below the stronger 
sources (plus symbols) for source sizes larger than ~ 100". This 
is because the insertion of the stronger sources in the astronom- 
ical source model used for iterations 1 to 20 helps to recover 
more astronomical signal when their outer parts start to reach the 
thresholds mentioned in Sect. 13.11 The spatial filtering is there- 
fore less severe for strong sources than for faint sources. These 
results are summarized in Table lA.ll 

A further Monte Carlo test was performed with elliptical 
Gaussian sources with a fixed aspect ratio of 2.5. The procedure 
to generate these artificial sources was the same, with the orien- 
tation of the sources also assigned randomly. The peak flux den- 
sities ranged between 85 and 500 mJy/beam and the size along 
the minor axis (FWHM) between 19.2" and 3.7'. Three sets of 
sources were analysed, providing a sample of 68 sources in to- 
tal. The results are very similar to the circular case, but with 
characteristic scales of the filtering along the minor axis of the 
elliptical sources 20-25% smaller (see Table lA.ll . In return, the 
corresponding characteristic scales of the filtering along the ma- 
jor axis are a factor of 2 larger than in the circular case. The 
peak flux densities are reduced by ~ 15% to ~ 20% for input 
sizes along the minor axis varying from ~ 180" to ~ 220" (i.e. 
sizes from ~ 450" to ~ 550" along the major axis). The aspect 
ratio is well preserved within ~ 10% after the last iteration, but 
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Fig. A.l. Statistical properties of 60 artificial, circular, Gaussian sources inserted into the raw time signals before data reduction (a 
to e) or directly inserted into the final unsmoothed continuum emission map of Cha I (f to j). a and f : ratio of output to input peak 
flux density as a function of input peak flux density, b and g: ratio of output to input peak flux density as a function of input size 
(FWHM). c and h: ratio of output major size to input size as a function of input size, d and i: ratio of output minor size to input 
size as a function of input size, e and j: ratio of output to input total flux as a function of input size. The size of the black symbols 
scales with the input size of the sources. In addition, crosses and plus symbols are for input peak flux densities below and above 150 
mJy/19.2"-beam, respectively. The red circles in panels b to e and g to j show a fit to the data points using the arbitrary 3-parameter 
function y = log(Qr/(l + (x//3y)). 



the sizes of the sources are reduced by ~ 15% to ~ 25% for in- 
put sizes along the minor axis from ~ 160" to ~ 220" (i.e. from 
~ 400" to ~ 550" along the major axis). Here again, sources 
with a peak flux density higher than ~ 150 mJy/beam are better 
recovered than weaker sources that cannot or can only partly be 
incorporated into the astronomical source model used at itera- 
tions 1 to 20. 

These two Monte Carlo tests clearly show that the reduction 
of LABOCA data based on correlations between pixels to reduce 
the atmospheric and instrumental noise can have a strong impact 
on the structure of the detected astronomical sources. Depending 
on their shape and strength, structures with (major) sizes larger 
than ~ 140" - 280" if circular or ~ 300" - 550" if elliptical 
with an aspect ratio of 2.5 are significantly filtered out (> 15% 
losses of peak flux density and size). Elongated structures are 
better recovered than circular ones. We note that the reduc- 
tion performed for data obtained at the Caltech Submillimeter 
Observatory (CSO) with the bolometer camera Bolocam suf- 
fers from the same problem, with a similar characteristic scale 
above which losses due to spatial filtering start to be significant 
(~ 120" for circular sources, see Enoch et al,:,2006) . 

Several attemps to attenuate the correlated noise removal by 
using the outer pixels only, or to suppress the correlated noise 
removal and replace it subscan-wise with a simple first-order 
baseline removal in the time signals of the OTF scans, were 
not successful. These alternative reduction schemes produced 
somewhat more extended emission in places where we expected 
some, but at the same time, there were much more atmospheric 
and instrumental noise residuals distributed within the map that 
also mimicked extended structures where we did not expect any, 
by comparison to the extinction and C'^O 1-0 maps for instance. 
Attempts to reduce the data without flattening the power spec- 




5 10 15 20 5 10 15 



Iteration Itei-tition 

Fig.A.2. a Convergence of the total 870 /vm flux recovered in 
Cha I by the iterative process of data reduction, b Relative flux 
variation between iterations / and ; - 1 as a function of iteration 
number i. 



trum at low frequency or with a lower cutofl" frequency did not 
bring any significant improvement either. 

A.2. Convergence of the iterative data reduction 

The iterative process of data reduction described in Sect. 13.1 1 al- 
lows to recover additional astronomical flux after each iteration. 
In this section, we study the convergence of this method. The 
evolution of the total 870 fim flux recovered in Cha I as a func- 
tion of iteration is shown in Fig. lA.2b . The relative flux variation 
is only 0.7% at iteration 20, which suggests that the data reduc- 
tion process indeed converged. 

However, the previous estimate is based on relative varia- 
tions of the astronomical flux from one iteration to the next one. 
It cannot completely exclude a very slow convergence. A more 
secure estimate can be derived from the set of artificial, circu- 
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Fig. A.3. Convergence of the iterative data reduction for artificial, circular, Gaussian sources with a peak flux density larger than 150 
mJy beam '. The first and second columns show the ratio of the output to the input peak flux densities, and its relative variations, 
respectively. The third and fourth columns show the same for the total flux. The different rows show sources of different widths, as 
specified in each panel of the first column. The size of the symbols increases with the input peak flux density. 



lar, Gaussian sources described in Appendix lA.ll Figure IA.3I 
demonstrates that the number of iterations needed to reach con- 
vergence depends on the size of the sources. We limit the anal- 
ysis to sources stronger than 150 mJy beam"' (fainter sources 
show the same behaviour as far as the relative variations from 
one iteration to the next one are concerned, but the curves are 
noisier). Convergence on the peak flux density and total flux with 
remaining relative variations < 2% is reached after a few itera- 
tions (3 to 5) for sources with FWHM < 160" (Figs. IA.3b . d, 
f, and h). Sources with 160" < FWHM < 320" reach conver- 
gence with relative variations < 2% at iteration 15 (Figs. IA.3I 
and 1). At iteration 20, the peak flux densities of sources with 
320" < FWHM < 380" have converged within 2% while the 
total fluxes may still vary by 3% (see Figs. IA.3h and p). The 
convergence for the single source with 380" < FWHM < 440" 
is good within 5% at iteration 20. The results are similar for the 
set of elliptical sources: convergence is reached within 2% at it- 
eration 4-10 for sources with size along the minor axis < 200", 
and within 2% at iteration 18 for sources with size ~ 220". 



This experiment shows that a significant number of iterations 
(about 20) are necessary to recover as much extended emission 
as possible. Enoch et al. (2006) have a similar conclusion with 
Bolocam for circular sources with FWHM < 120" (conver- 
gence at iteration 5, see their Fig. 4), but they do not mention 



or show how many iterations they need to reach convergenc^B 
for much larger sources. 

Appendix B: Assumptions 

In this appendix, we explain the various assumptions that are 
used to analyse and interpret the Cha I data. We also explain 
how the column densities and average densities obtained in the 
past for other molecular clouds are rescaled in the present study 
to compare with Cha I (Table |7]). 

B.1. Distance 

IWhittet et al.l d 19971) analysed the distribution of extinctions of 
field stars with distance along the line of sight of Cha I and de- 
rived a distance of 150 + 15 pc for this cloud. They also de- 
duced a photometric distance for the B-type star HD 97300 of 
152 ± 18 pc. They combined these results with the parallaxes of 
this and another star (HD 97048) measured with the Hipparcos 
satellite with a larger uncertainty (190 ± 40 and 180 + 20 pc, 
res pectively) and adopted a "final" distance of 160 ± 15 pc. 

iKnude & H^gl (Il998h used the interstellar reddening from 
the Hipparcos and Tycho catalogues to derive the distance to 



Convergence for artificial sources means here to reach a stable flux 
density, not to reach a flux density equal to the input one, which is never 
the case for large sources that will always suffer from spatial filtering 
due to the correlated noise removal. 
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Cha I. They found a distance of 150 pc as revealed by a sharp 
increase in their plot of E(B-V) color excess versus distance. 

However, Knude (2 010) recently derived a further distance 
of 196 + 13 pc for Cha I from a fit to an extinction vs. distance 
diagram. The extinction was obtained from a (H-K) vs. (J-H) 
diagram from 2MASS and the distances were estimated from a 
(J-H)() vs. Mj relation based on Hipparcos. Rather than consid- 
ering the distance at which the first stars show a steep increase 
in extinction as the distance of the cloud, they fitted a predefined 
function to a subsample of the data. It is unclear to us whether 
this fitting procedure does give a more reliable distance estimate 
because the selection of the subsample defining the extinction 
jump for the fit may introduce a bias. Therefore, as long as this 
further distance has not been confirmed by another method, we 
keep the older estimate (150 pc) as the distance of Cha 1. If the 
further distance turns out to be correct, then all physical sizes and 
masses derived in this article are underestimated by a factor 1.3 
and 1.7, respectively, and the mean densities are overestimated 
by a factor 1.3. The column densities are not affected. 

B.2. Temperature 



iToth et al.1 (l2000b derived a dust temperature of 12 + 1 K for 
the "very cold cores" they identified in Cha 1 on somewhat 
large scales (HPBW ~ 2') from their analysis of the ISOPHOT 
Serendipity Survey (see cores 1 to 6 in their Table 2). The dust 
temperature is most likely not uniform and this value rather rep- 
resents a mean temperature. The dust temperature is expected 
to decrease toward the center of starless dense cores that do 
not contain any central heating, even down to ~ 5.5 - 7 K as 
measured for the gas in the prestellar cores L1544 and L183 by 
[Crapsi et al. (2007) and Pagani et al. (2007), respectively. The 
modeled dust temperature profil es of the pre s tellar cores L1512, 
L1689B, and L1544 derived bv Evans et al.1 (1200 ih indeed vary 
from 12- 13 K in the outer parts to 7-8 K at the core center. 
Similarlv. lNutter et al.l (12008 ^ measured a dust temperature vary- 
ing from 12 K in the outer parts to 8 K in the inner parts of a 
filament in the Taurus molecular cloud 1 (TMCl). Thanks to its 
good angular resolution and its good sampling of the spectral 
energy distribution of these cores, such temperature variations 
will be routinely measured in the far-infrared with the Herschel 
satellite, which wil l help refining the mass a nd column density 
estimates (see, e.g.. lAndre & Saracenol2005h . For simplicity, we 
assume here a uniform dust temperature of 12 K. 

If the temperature inside the dense cores drops to, e.g., 7 K, 
their masses and column densities computed in this article could 
be underestimated by up to a factor of 3, depending on the exact 
shape of the temperature profile. However, the total mass of a 
core being dominated by its outer parts, it should not be that 
much underestimated. 

B.3. Conversion of visual extinction to H2 coiumn density 

The standard factor to convert a visual extinction Ay to an H2 
column density is 9.4 x 10^° cm"^ mag"'. It is based on the ra- 
tio Nh/E(B - V) ^ 5.8 X 10^' atoms cm ^ mag ' measured by 
iBohlin et alj (11978!) with data from the Copernicus satellite, and 
a ratio of total-to-selective extinction Ry = Av/E{B - V) - 3.1 
corresponding to the standard extinction law that is represen- 
tative of dust in diffuse regions in the local Galaxy (see, e.g., 
[Oraine 2003 and references therein). 

Several measurements suggest tha t Ry is larger toward 
denser molecular clouds. In particular, IChapman et al.l (l2009l) 



computed the mid-lR extinction law based on Spitzer data in 
the three molecular clouds Ophiuchus, Perseus, and Serpens 
and found indications of an increase of Ry from 3.1 for Ak^ < 
0.5 mag to 5.5 for A^^ > 1 mag, which they interpret as evidence 
for grai n growth. U sing the parameterized extinction law of 
Cardelli et al.1 (Il989l their equations 1 , 2a, and 2b) with Ry ^ 3.1 
and 5.5, respectively, we find that these two thresholds corre- 
spond to Ay < 4.2 mag and Ay > 7.2 mag. Vrba & Rvdgre^ 
(il984 came to a similar conclusion for Cha I since they found 
a standard extinction law toward the outer parts of the cloud 
but reported a larger Ry ( > 4) for field stars near the denser 
parts of the cloud (see also Whi ttet et alj 19941). Theref o re, we 
assume that the thresholds derived by Chapman et alJ (l2009l) 
also apply to Cha 1. With Ry = 5.5, the factor to convert 
a visual extinction in to an H2 column density becomes 6.9 X 
10^" cm"^ mag ' (see lE vans et al.ll2b09l and references therein). 
However, the visual extinction values presented here are de- 
rived from near-infrared measurements assuming the standard 
extinction law {Ry - 3.1). Compensating for this assumption, 
the "corrected" conversion factor for Ry - 5.5 mag becomes 
Nuj/Ay = 5.8 X 10^" cm"^ mag"' for our extinction map cali- 
brated with Ry - 3.1. 

In summary, for Ay below ~ 6 mag in our extinction map we 
assume Ry - 3.1 and the standard conversion factor N-n^jAy = 

9.4 X 10^*^ cm"^ mag"'. For Av above ~ 6 mag, we assume Ry = 

5.5 and Nh^/ Ay = 5.8 x 10^" cm^^ mag"'. 



B.4. Dust mass opacity 

The dust mass opacity in the submm range is not well con- 
strained and is typically u ncertain by a factor of 2 (e.g. 
l Ossenkopf & HenningI 119941) . Since most sources detected in 
the survey of Cha I are weak, starless cores, we adopt a dust 
mass opacity, /cgvo, of 0.01 cm- g"' in most cases (and a dust 
opacity in dex /3 = 1.85) . It co rresponds to the value recom- 
mended bv lHenning et al.l (Il995l) for "cloud envelopes", assum- 
ing /3 ~ 1.5 - 2.0 and a gas-to-dust mass ratio of 100. For higher 
densities (ne^ ~ 5 x 10^ cm"-^) that characterize the inner part 
of a few sources in ou r survey. Km) is most likely h igher, on the 
order of 0.02 cm^ g"' (lOssenkopf & Henningll994l) . Finally, the 
dust mass opacity in c ircu mstellar disks is expected to be even 
higher. We follow Beckwith et al.l d 19901) with (3=1 and derive 
^■870 = 0.03 cm^ g"' for this type of sources. In the article, we 
each time specify which dust mass opacity we use to compute 
column densities and masses. We assume in all calculations that 
the dust emission is optically thin. 

B.5. Moiecuiar weight 

The mean molecular weight per hydrogen molecule and per free 
particle, yUe, and yu, are assumed to be 2.8 and 2.37, respectively 
(see, e.g., KaufFm ann et al. 2008). The former is used to con- 
vert the flux densities to H2 column densities while the latter is 
needed to estimate the free-particle densities. 

B.6. Rescaiing of previous results on other clouds 
B.6.1. Column densities 

Assuming a temperature of 10 K and a dust m ass opacity of 
KiA mm = 0.0114 cm^ g"'. lEnoch et al ] (I2008h derive average 
peak H2 column densities of 1.2 and 1.0 xlO^^ cm"^ in Perseus 
and Serpens, respectively. Their values have to be multiplied by 
1.3 to match the same properties as we use for Cha 1 (12 K, 



33 



A. Belloche et al.: The end of star formation in Cliamaeleon I ? 



j6 = 1.85 yielding a-268GHz = 0.0063 cm^ g"'). In addition, 
both clouds are at larger distances than Cha I and the Bolocam 
maps have a poorer angular resolution (31") compared to our 
LABOCA map (21.2"). If the sources in Perseus and Serpens 
are not uniform in the Bolocam beam but rather centrally peaked 
with a column density scaling proportionally to the inverse of 
the beam size, then their peak column densities measured in the 
same physical size as for Cha I would be about 2.4-2.5 times 
larger The full rescaling factor is thus 3.1-3.5. We rescale the 
av erage colum n density of the SCUBA starless cores analysed 
bv lCurtis & R icher (2010) in Perseus in the same way (12 K and 
/("sssGHz = 0.010 cm- g"'), after correction for a small errof^ in 
their numerical equation 5. The column density is also rescaled 
proportionally to the inverse of the beam size. The rescaled val- 
ues of both studies of Perseus are consistent. 

The average peak column densi ty of the Bolocam starless 
cores observed by Enoch et al] ((2008 ) in Ophiuchus was rescaled 
in the same way as above. The average peak colum n density de- 
rived for the starless cores of lMotteet aP (Il998h after rescal- 
ing to the same physical size as for the Cha I sources with an 
Ophiuchus distance of 125 pc is a factor of ~ 2 lower The dis- 
crepancy between both studies may come from the background 
emission not being removed in the former study, from the as- 
sumption of a unique temperature for all sources in the former 
study while the second one uses different temperatures (from 12 
to 20 K), and from our assumption of a column density scaling 
wit h the inverse of the radiu s . 

iKauffmann et alJ (l2008l) analysed their maps of Taurus after 
smoothing to 20", i.e. nearly the same spatial resolution as our 
Cha I map given also the slightly different distances. Their col- 
umn density sensitivity is similar to ours. Assuming the same 
dust properties as for Cha I (T^usi = 12 K and = 1.85 yield- 
ing /^asoGHz = 0.0055 cm^ g^), the average column density in 
an aperture of 20" derived from the flux densities listed in their 
Table 5 is 2.3 j x 10^^ cm"^ for their sample of 28 Taurus star- 
less "peaks". 

B.6.2. Densities 

lEnoch et"an (12008) obtain average densities of 1.7, 1.4, and 2.8 x 
10^ cm"^ within an aperture of diameter 10"* AU for the starless 
cores in Perseus, Serpens, and Ophiuchus, respectively. We as- 
sume a density profile proportional to r with p ~ 1 .5-2, which 
is typical for such cores, and rescale their densities with a tem- 
perature of 12 K and a dust opacity /caasGHz = 0.0063 cm^ g"'. 

With th e same assumptions as in Sect. IB. 6. II the 28 starless 
"peaks" of IKauffmann et al] ( l2008l) in Taurus have an average 
mass of 0.51^Q3 within an aperture of diameter 8400 AU. 
We assume the same density profile as above to compute the 
average mass within an aperture of diameter 7500 AU and derive 
the average density. 



Appendix C: Multiresolution decomposition 

As mentioned in Sect. 1.5.1 of IStarck et al.l (Il998h . multireso- 
lution tools based on a wavelet transform suffer from several 
problems, in particular they create negative artefacts around the 
sources. In addition, since the emission is smoothed with the 
scaling function (O) to produce the large-scale planes, the shape 
of the large-scale structures is not preserved and tends toward 



Table C.l. Characteristics of the multiresolution decomposi- 
tion as measured on a sample of circular and elliptical artificial 
Gaussian sources with a signal-to-noise ratio higher than 9. 



Input size 




Scale at which 






/;//; > 40% 


fp/f;, > 80% 


~ 100% 


(1) 


(2) 


(3) 


(4) 


< 30" 


2 


3 


4 


30-60" 


3 


4 


5 


60-120" 


4 


5 


6 


120-200" 


5 


6 


7 


200-300" 


6 


7 


7 


300-400" 


6 


7 





Notes. The input size is the FWHM for circular sources, and 
^JFWHMI X FWHM2 for elliptical sources, with FWHMl and 
FWHM2 the sizes along the major and minor axes, respectively. f° 
and fp are the output and input peak flux densities, respectively. 



the shape of O (which, in turn, creates the negative artefacts in 
the small-scale planes mentioned above). For instance, a circu- 
larly symmetric filter will tend to produce roundish structures on 
large scale. To avoid these issues, we performed a multiresolu- 
tion anal ysis based on the m edian transform following the algo- 
rithm of Stocketal] ([1998', Sect. 1.5.1) that we implemented in 
Python. In short, the input image is first "smoothed" with a fil- 
ter that consists of taking the median value on a 3 x 3-pixel grid 
centered on each pixel. The difference between the input and 
smoothed images yields the "detail" map at scale 1 (the smallest 
scale). The same process is repeated with a 5 x 5-pixel grid using 
the smoothed map as input map, and the difference between the 
two images produces the detail map at scale 2. The process is 
done 7 times in total, yielding detail maps at 7 different scales. 
The size of the convolution grid at scale / is « x n pixels with 
n = 2' -H 1, i.e. 18.3", 30.5", 54.9", 1.7', 3.4', 6.6', and 13.1' for 
scales 1 to 7, respectively. The last smoothed map contains the 
residuals (very large scale). The sum of the residual map and the 
7 detail maps is strictly identical to the input map. 

The properties of this multiresolution decomposition were 
estimated with the samples of artificial Gaussian sources of 
Appendix |AT| (both circular and elliptical with aspect ratio 2.5) 
inserted into the reduced continuum map of Cha I. Table IC.ll 
presents the results for sources with a signal-to-noise ratio 
higher than 9. For each scale i, the artificial sources were fit- 
ted in the map containing the sum of the detail maps from 
scale 1 to / (hereafter called sum map at scale /). Compact 
sources (FWHM < 30") are fully recovered at scale 4 while 
sources with FWHM > 120" are strongly filtered out (losses 
> 60%). Therefore, scale 4 is a good scale to extract com- 
pact sources when they are embedded in larger scale cores 
(> 120"). Similarly, scale 5 is a good scale to extract sources 
with FWHM < 120" with an accuracy better than 20% on 
the peak flux density, while significantly filtering out the more 
extended background emission (losses > 60% for FWHM > 
200"). In the case of the elliptical sources with aspect ratio 
2.5, the conclusions are the same, the key parameter being the 
geometrical mean of the sizes along the major and minor axes 
( FWHMl X FWHMl). 



The solid angle of a Gaussian beam is = 4 \Ip) ^hpbw^ ^i\h Ohpbw 
the ha l f-power beam width (se e, e.g.. Appendix A.3 of lKauffmann et afl 
l2008h . lCurtis & Richeil (l2010h seem to use = j&lfpg-^- 



List of Objects 

'Cha-MMS2' on pageE 
'Cha-MMSr on page 4 
'WW Cha' on page^"^ 
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'IRAM 04191 + 1522' on page|25| 
'L1521F' on pagejH 
'L1448 IRS2E' on pageHH 



